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I.  INTRODUCTION 


Formaldehyde  is  a  ubiquitous  flame  zone  species.  Reactions  between 
formaldehyde  and  various  oxidizers  provide  a  major  portion  of  the  heat 
released  in  the  primary  flame  zone  of  nitrate  esters  and  nitramines.l 
Any  model  or  description  of  propellant  combustion  must  account  for  the 
oxidation  of  formaldehyde.  A  precise  model  requires  detailed  modeling 
of  the  flame  chemistry,  and  thus  rate  coefficients  for  chemical  reactions. 
Since  experiment  cannot  measure  rate  coefficients  for  all  the  elementary 
flame  reactions,  it  is  important  to  develop  and  test  methods  for  esti¬ 
mating  those  rate  coefficients. 

Thermochemical  kinetics,  described  and  justified  in  Benson's  mono¬ 
graph,  provides  one  approach  to  estimate  rate  coefficients  for  several 
classes  of  chemical  reactions. 2  These  methods  require  thermochemical 
parameters,  such  as  heat-of-formation,  entropy  and  heat  capacity,  for 
each  chemical  species  involved  in  the  chemical  reaction.  For  many  com¬ 
pounds,  these  data  can  be  found  in  the  various  compendia  of  thermochemical 
data;  but  for  other  species,  accurate  thermochemical  data  are  not  avail¬ 
able.  Furthermore,  to  estimate  rate  coefficients  for  unimolecular  re¬ 
actions  of  radical  species  requires  thermochemical  data  for  evanescent 
transition-state  species;  data  not  susceptible  to  direct  measurement. 

This  report  describes  the  application  of  many-body  perturbation  theory 
to  predict  thermochemical  parameters  for  some  flame  zone  reactions  that 
involve  formaldehyde,  and  for  several  other  species  that  are  important 
in  models  that  include  formaldehyde. 

The  chemical  reactions  to  be  examined  are: 


CH3OH  -►  CH30  +  H  (I) 

CH30  CH20  +  H  (II) 

CHo0  -»■  CHO  +  H  (HI) 

CHO  -»•  CO  +  H  (IV) 

CH20  -v  H2  +  CO  .  (V) 

CH20  HCOH.  (VI) 


1 

R.  A.  Fifer^  "High  Temperature  Pyrolysis  of  Methyl  and  Ethyl  Nitrate "3 
Seventeenth  Combustion  Institute  Symposium ,  University  of  Leedss  England , 
1978. 


2 


S. 


W.  Benson ,  Kinetics ,  John  Wileyi  New  York 3  1976. 
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This  selection  implies  neither  that  all  of  these  reactions  are  important 
flame-zone  reactions  nor  that  other  reactions  are  unimportant.  The 
selected  reactions  represent  only  a  fraction  of  the  many  flame-zone  re¬ 
actions  that  involve  formaldehyde. 

Reactions  (I)  and  (II)  both  include  the  methoxy  radical,  CH3O.  Ac¬ 
curate  predictions  of  the  heat-of-reaction  for  each  of  these  processes 
should  provide  enough  information  to  determine  the  heat-of-formation  for 
this  radical.  Several  unimolecular  reactions  that  remove  formaldehyde 
from  a  flame  are  given  by  equations  (III),  (V)  and  (VI),  while  reaction 
(IV)  competes  with  the  bimolecular  reaction  between  the  formyl  radical 
and  formaldehyde.  For  each  of  the  listed  reactions,  theoretical  predic¬ 
tions  of  the  heat-of-reaction  are  obtained.  Activation  energy  barriers 
are  estimated  for  reactions  (IV),  (V)  and  (VI). 

The  theoretical  methods  employed  in  this  study  are  described  in  the 
next  section  and  in  several  appendices.  Succeeding  sections  describe 
details  of  the  electronic  structure  calculations  and  results  of  the 
thermochemical  estimates  for  each  reaction;  the  final  section  discusses 
the  implications  for  chemical  models  of  the  flame  zone. 


II.  DISCUSSION  OF  THEORETICAL  METHODS 

For  many  years  electronic  structure  theorists  have  attempted  to 
derive  simple  schemes  for  predicting  or  calculating  thermochemical  para¬ 
meters  for  chemical  species.  The  best  known  of  the  models,  MINDO,  de¬ 
veloped  by  Dewar  and  co-workers,  predicts  thermochemical  data  with  fair 
accuracy  for  hydrocarbons.3  However,  the  performance  of  the  method 
diminishes  for  non-hydrocarbon  organic  species  because  there  are  insuf¬ 
ficient  empirical  data  to  develop  fitting  parameters.  In  practice,  this 
means  that  the  method  fares  poorly  when  used  to  study  oxygen  containing 
compounds . 

When  Snyder  and  Basch  used  ab  initio  molecular  orbital  methods  to 
study  the  heats-of-reaction  of  closed-shell  molecules,4  they  obtained 
fair  agreement  between  theory  and  experiment.  Although  improved  agree¬ 
ment  might  result  from  an  increase  in  the  size  of  the  atomic  basis  sets 
used  in  the  calculations,  recent  research  suggests  that  more  fundamental 
problems  must  be  addressed.  In  particular,  the  accurate  prediction  of 
relative  energies  requires  accounting  for  the  correlation  between 


g 

M.  J.  S.  Dewar  and  R.  C.  Doughety  3  "The  PMO  Theory  of  Organic  Chemistry" 3 
Plenum  Press 3  New  York 3  1975. 

4 

L.  C.  Snyder  and  H.  Basch3  "Heats  of  Reaction  from  Self-Consistent  Field 
Energies  of  Closed-Shell  Molecules "3  J.  Am.  Chem.  Soc.  91 3  2189  (1969). 
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electron  motions,  an  effect  that  is  not  treated  in  self-consistent  field 
calculations .5,6,7 

Historically,  estimates  of  molecular  correlation  energy  were  obtain¬ 
ed  by  using  the  method  of  configuration  interaction  (Cl)*  to  improve  the 
wavefunction  obtained  from  a  molecular  orbital  (MO)  SCF  calculation.  In 
the  past  several  years,  theoretical  techniques  based  upon  the  linked -dia¬ 
gram  theorem  have  been  extended  to  molecular  applications,  and  several 
comparisons  of  the  various  theoretical  approaches  to  the  correlation 
problem  have  been  published. 5, 8  jn  recent  works,  Pople  suggests  a  set 
of  desiderata  for  computational  techniques: 

(a)  The  method  should  be  size-consistent.  An  energy  calculated 

for  an  assembly  of  isolated  molecules  should  equal  the  sum  of  the  energies 
calculated  for  individual  molecules. 

(b)  The  method  should  be  invariant  under  transformations  within  a 
set  of  degenerate  MO's. 

(c)  The  method  should  be  variational.  The  energy  should  be  an 
upper-bound  to  the  energy  obtained  by  an  exact  solution  to  the  electronic 
Schrodinger  equation. 

Few  theoretical  methods  satisfy  all  three  criteria  simultaneously.  Spin- 
unrestricted  Hartree-Fock  theory  (UHF)  often  satisfies  all  three  criteria, 
but  the  method  takes  no  account  of  correlation  effects  between  electrons 
of  different  spins.  A  complete  solution  of  the  configuration  interaction 
problem  for  a  given  basis  set  satisfies  all  three  criteria,  but  such  cal¬ 
culations  are  practical  only  for  small  molecular  systems.  As  pointed  out 
by  Pople,  there  is  a  need  for  theories  at  an  intermediate  level  that 
account  for  electron  correlation. 5 


J.  A.  Pople ,  R.  Krishnan ,  H.  B.  Sohlegel  and  J.  S.  Binkley ,  "Electron 
Correlation  Theories  and  Their  Application  to  the  Study  of  Simple  Reaction 
Potential  Surfaces ",  Int.  J.  Quantum  Chem.  14,  545  (1978). 

g 

R.  J.  Bartlett,  I.  Shavitt  and  G.  D.  Purvis,  "The  Quadratic  Force  Field 
of  RgO  Determined  by  Many-Body  Methods  that  Include  Quadruple  Excita¬ 
tion  Effects",  submitted  to  J.  Chem.  Phys . 

7 

J.  D.  Goddard  and  H.  F.  Schaefer,  III,  "The  Photodissociation  of  For¬ 
maldehyde:  Potential  Energy  Surface  Features ",  J.  Chem.  Phys.  70,  5117 
(1979).  y 

g 

R.  J.  Bartlett  and  G.  D.  Purvis ,  "Many-Body  Perturbation  Theory ,  Coupled- 
Pair  Many -Electron  Theory,  and  the  Importance  of  Quadruple  Excitations 
for  the  Correlation  Problem",  Int.  J.  of  Quantum  Chem.  14,  561  (1978). 

* 

A  description  of  the  method  of  configuration  interaction  is  given  in 
Appendix  A. 
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The  technique  most  often  used  to  estimate  correlation  energy  is  a 
truncated  Cl  calculation.  Generally,  the  many  electron  wave  function 
consists  of  the  Hartree-Fock  configuration,  $0,  plus  all  configurations 
obtained  by  double-excitations  of  electrons  from  one-electron  functions 
occupied  in  <J>0  to  the  virtual  orbitals  of  $o-*  A  configuration  inter¬ 
action  calculation  that  uses  such  a  many-electron  wave  function  is  called 
configuration-interaction  with  double-excitations  (CID) .  CID  calculations 
satisfy  criteria  (2)  and  (3) ,  but  do  not  satisfy  the  size-consistency 
requirements,  since  the  wave-function  cannot  account  for  simultaneous 
double-excitations  in  different  molecules. 


Among  alternatives  to  conventional  Cl  are  theoretical  techniques 
based  on  the  linked-diagram  theorem. 9.10  One  of  these  techniques  is  the 
couple-cluster  approach  (CCA) . ^ ^  When  the  cluster  analysis  is  limited 
to  a  consideration  of  only  two-body  terms,  which  include  both  the  double¬ 
excitations  included  in  CID  and  the  simultaneous  double-excitation  omitted 
by  the  CID,  the  method  is  called  coupled  cluster  with  double  substitutions 
(CCD),5  or  coupled-pair  many-electron  theory  (CPMET).13  Relative  to  a 
Hartree-Fock  reference  function,  this  theory  satisfies  criteria  (1)  and 
(2),  but  does  not  satisfy  the  variational  requirement.** 


Another  linked-diagram  related  approach  is  many-body  perturbation 
theory  (MBPT).5  This  technique  uses  the  Moller-Plesset  partitioning  of 


g 

K.  A.  Brueckner,  "Two-Body  Forces  and  Nuclear  Saturation.  III.  Details 
of  the  Structure  of  the  Nucleus" ,  Phys.  Rev.  97,  1353  (1955). 


10 

K.  A.  Brueckner,  "Many-Body  Problem  for  Strongly  Interacting  Particles. 
II.  Linked  Cluster  Expansion ",  Phys.  Rev.  100,  36  (1955). 


'J.  Cizek,'  "On  the  Correlation  Problem  in  Atomic  and  Molecular  Systems. 
Calculation  of  Wavefunction  Components  and  Ursell-Type  Expansion  Using 
Quantum  Field-Theoretical  Methods ",  J.  Chem .  Phys.  45,  4256  (1966). 

J.  PalduSj  J .  Cizek  and  I.  Shavitt 3  " Correlation  Problems  in  Atomic 
and  Molecular  Systems .  IV.  Extended  Coupled  Pair  Many-Electron  Theory 
and  its  Application  to  the  BH$  Molecule" 3  Phys.  Rev.  A5>  50  (1972). 

? 

J .  Cizek,  "On  the  Use  of  the  Cluster  Expansion  and  the  Technique  of 
Diagrams  in  Calculations  of  Correlation  Effects  in  Atoms  and  Molecules" , 
Adv.  Chem.  Phys.  14,  35  (1969). 


Virtual  orbitals  are  those  one-electron  functions  that  are  not  occupied 
in  the  reference  function  $Q. 

t 

A  description  of  coupled-pair  many-electron  theory  is  given  in  Appendix 

B. 
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the  molecular  Hamiltonian,^  such  that  the  Hartree-Fock  solution  is  taken 
to  be  the  unperturbed  Hamiltonian,  H0,  and  the  residual  part  of  the  mole¬ 
cular  Hamiltonian  is  treated  as  the  perturbation,  V.  Bartlett  and  co¬ 
workers  have  pioneered  molecular  applications  of  many-body  perturbation 
theory, 8,15  a  technique  that  satisfies  desiderata  (1)  and  (2),  but  like 
CCD,  does  not  satisfy  a  variational  principle.* 

Reliable  predictions  of  potential  energy  hypersurfaces  and  thermo¬ 
chemical  parameters  require  that  computational  methods  exceed  the  Hartree- 
Fock  self-consistent-field  level  of  accuracy  and,  thus,  must  include  cor¬ 
relation.  The  choice  of  method  used  to  estimate  the  correlation  energy 
in  such  studies  depends  upon  the  importance  of  size-consistency  effects 
and  the  related  size-extensivity  criterion.**  Whenever  these  effects 
are  important,  truncated  Cl  methods,  such  as  CID,  do  not  give  accurate 
representations  of  potential  energy  hypersurfaces.  Size-extensive  meth¬ 
ods,  such  as  CCD  and  MBPT,  provide  more  uniform  representations  of  po¬ 
tential  energy  hypersurfaces  than  does  a  truncated  Cl  calculation. 

Bartlett,  et  al . ,  concluded  that  the  effect  of  size-consistency  on  elec¬ 
tronic  structure  calculations  was  important,  even  for  calculating  the 
potential  energy  near  the  equilibrium  structure  for  a  molecule  as  small 
as  water. 6  Since  size-consistency  and  size-extensivity  depend  upon  the 
number  of  electrons,  theoretical  methods  that  satisfy  these  criteria 
should  predict  more  accurate  values  for  thermochemical  parameters  than  a 
method  such  as  CID.  Consequently,  this  study,  in  the  main,  uses  fourth- 
order  MBPT  to  estimate  relative  energies.  The  MBPT  calculations  are 
supplemented  by  CCD  calculations  for  some  crucial  molecular  configurations. 
The  latter  calculations  serve  to  check  the  convergence  of  the  MBPT  calcu¬ 
lations  .8 

Several  appendices  discuss,  in  more  formal  terms,  the  theoretical 
methods  mentioned  above,  and  provide  a  discussion  of  some  aspects  of 
size-consistency  and  size-extensivity  in  molecular  calculations.  Most 
of  the  information  presented  in  the  appendices  can  be  found  in  various 
papers  by  Bartlett,  et  al. 6,8,15 

III.  DETAILS  OF  THE  CALCULATIONS 

The  electronic  states  of  the  molecules  discussed  in  this  report  are 
represented  at  the  Hartree-Fock  level  by  unrestricted  Hartree-Fock  (UHF) 


14  C.  Moller  and  M.  S.  Plesset 3  "Note  on  an  Approximation  Treatment  for 
Many-Electron  Systems* r,  Phys.  Rev .  46>  618  (1934). 

15 R.  J.  Bartlett  and  D.  M.  Silver, 3  nSome  Aspects  of  Diagrammatic  Per¬ 
turbation  Theory” >  Int .  J.  of  Quantum  Chem.  9S 3  183  (1975). 

*A  description  of  many-body  perturbation  theory  is  given  in  Appendix  C . 

*  *  .  . 

Size-consistency  and  size-extensivity  are  discussed  Appendix  D. 
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wavefunctions .  For  closed  shell  molecules,  such  as  formaldehyde  and 
carbon  monoxide,  these  wavefunctions  are  just  those  that  would  be  obtain¬ 
ed  with  a  restricted  Hartree-Fock  (RHF)  calculation.  As  noted  in  Section 
II,  a  UHF  calculation  usually  satisfies  all  of  Pople’s  desiderata. 

The  computations  reported  here  refer  to  a  double-zeta  plus  polariza¬ 
tion  atomic  basis  set  consisting  of  contracted  Gaussian-type  functions.* 
Dunning’s  4s3pl^  contraction  of  Huzinaga’s  9s5p  primitive  setl?  was  used 
for  the  first-row  atoms.  The  carbon  and  oxygen  basis  sets  were  each  aug¬ 
mented  by  a  single  set  of  d-orbital  polarization  functions  with  exponents 
0.75  and  0.85,  respectively.  The  hydrogen  basis  set  included  the  3s  con¬ 
traction  of  Huzinaga’s  4s  primitive  set,  with  scaled  exponents  (6  =  1.2), 
augmented  by  a  p-type  polarization  function  with  orbital  exponent  1.0. 

All  integrals  over  atomic  basis  functions  were  calculated  with  the  MOLE¬ 
CULE  integral  program. 18 

The  SCF,  MBPT  and  CCD  calculations  were  performed  using  the  program 
PRPGTR. 18  The  different  calculations  include  SCF;  fourth-order  MBPT 
including  all  double-excitation  diagrams  that  contribute,  D-MBPT(4) ; 
fourth-order  MBPT  including  all  single-  and  double-excitation  diagrams, 
SD-MBPT (4) ;  fourth-order  MBPT  including  all  single-,  double-,  and  quad¬ 
ruple-excitation  diagrams  that  contribute,  SDQ-MBPT(4);  and  coupled-pair 
many-electron  theory,  CCD.  Although  relative  energies  reported  repre¬ 
sent  energy  differences  between  calculations  using  the  full  MO-basis, 
molecular  structure  parameters  were  optimized  at  SD-MBPT(4)  level  with 
the  restriction  that  the  core-orbitals  were  not  correlated  with  the 
remainder  of  the  MO-basis.  Molecular  structure  optimization  was  perform¬ 
ed  for  all  molecular  configurations  reported,  except  the  transition 
states  of  formaldehyde  reactions  leading  to  (a)  molecular  products,  and 
(b)  hydroxycarbene .  In  these  two  cases,  we  use  transition  state  geomet¬ 
ries  predicted  by  Goddard  and  Schaeffer. ? 


T .  H .  Dunning s  Jr .,  "Gaussian  Basis  Sets  for  Use  in  Molecular  Calcu¬ 
lations .  I.  Contraction  of  (9s >  Sp)  Atomic  Basis  Sets  for  First-Row 
Atoms ",  J .  Chem.  Phys.  53,  2823  (1970). 


S.  Huzinaga,  "Gaus sian-Type  Functions  for  Polyatomic  Systems. 
Chem.  Phys .  42,  1293  (1965). 


I. 


tt 

j 


J. 


1 8 

MOLECULE  was  written  by  J.  Almlof.  The  program  -is  described  in  "Pro¬ 
ceedings  of  the  Second  Seminar  on  Computational  Problems  in  Quantum 
Chemistry" 3  Strasbourg 3  France 3  1972  (1973).  PRPGTR  includes  program 
GRNFNCj  written  by  G.  D.  Purvis  and  program  UMBPT3  written  by  R.  J. 
Bartlett  and  G.  D.  Purvis. 


This  is  the  smallest  basis  set  that  can  provide  predictions  of  chemical 
accuracy . 
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IV.  DISSOCIATION  OF  METHANOL:  CH3OH  +  ClUO  +  H 

The  dissociation  of  methanol,  CH3OH,  to  radical  products  methoxy, 

CH3O,  and  hydrogen,  H,  is  not  generally  important  in  flame  zone  chemistry 
schemes  except  for  the  pyrolysis  of  methanol.  This  reaction  offers  a 
mechanism  to  estimate  the  heat-of-formation  of  the  methoxy  radical,  hence 
the  inclusion  in  this  study.  Previous  research  by  Redmon,  et  al.,^®  and 
Adams,  et  al.,^0  demonstrated  that  MBPT  predicts  heats-of-reaction  and 
dissociation  energies  that  agree  well  with  available  experimental  values. 
Since  heat-of-formation  data  exist  for  methanol  and  hydrogen,  the  heat- 
of-formation  of  methoxy  can  be  calculated  if  the  theory  predicts  an 
accurate  value  for  the  dissociation  energy  of  methanol.* 

To  predict  the  heat-of-reaction  for  the  dissociation  of  methanol,** 
electronic  energies  must  be  calculated  for  methanol,  methoxy  radical  and 
the  hydrogen  atom.  Since  linked-diagram  related  methods,  such  as  MBPT 
and  CCD,  are  size-consistent,  the  electronic  energies  of  the  products  are 
additive,  so  that  separate  calculations  are  performed  for  each  of  the 
products.  This  feature  offers  two  benefits  in  the  calculations  for 
methanol.  First,  two  separate  calculation  save  computer  time  and  reduce 
required  computer  resources  when  compared  to  the  supermolecule  calculation. 
Second,  separation  of  the  product  channel  to  two  calculations  permits  a 
thorough  analysis  of  the  Jahn-Teller  energy  lowering  for  the  methoxy 
radical . 

A  theoretical  study  of  methoxy  by  Yarkony,  et  al.,  showed  that  the 
ground  electronic  state  for  the  molecule  in  a  configuration  with  C3V  sym¬ 
metry  belonged  to  the  two-dimensional  (E)  irreducible  representation  of 
the  point  group. 21  The  Jahn-Teller  theorem  states  that  a  molecule  cannot 


19 

L •  T.  Redman 3  G.  D .  Purvis  and  R .  J.  Bartlett 3  nThe  Unimolecular 
Isomerization  of  Methyl  Isooyanide  to  Methyl  Cyanide” 3  J .  Chem .  Phys . 
69,  5386  (1978). 


20 

G .  F .  Adams 3  G .  D.  Bent3  G.  D.  Purvis  and  R.  J.  Bartlett 3  nThe  Elec¬ 
tronic  Structure  of  the  Formyl  Radical 3  HCO,r3  to  be  published 3  J.  Chem. 
Phys.3  1979 . 


21 

D.  R.  Yarkony 3  H.  F.  Schaefer  III.  and  S.  Rothenberg3  n Geometries  of 
the  Methoxy  Radical  (X2E  and  A2Aj  States)  and  the  Methoxide  Ionn3  J. 
Am.  Chem.  Soc.  96 656  (1974). 


o 

It's  clear  that  tH-^(CH^O)  can  also  be  estimated  by  predicting  the  dis¬ 
sociation  energy  of  the  radical  to  form  formaldehyde  and  hydrogen  and 
combining  this  result  with  known  heats-of- formation  of  the  products. 
This  calculation  will  be  described  in  the  next  section. 


The  heat-of-reaction  for  a  dissociation  reaction  is  often  called  the 
dissociation  energy.  Both  expressions  will  occur  in  this  report. 
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exist  in  a  degenerate  electronic  state,  requiring,  therefore,  a  distor¬ 
tion  of  the  molecular  configuration  to  destroy  the  spatial  symmetry  of 
the  state.  Yarkony,  et  al . ,  considered  a  single  vibrational  mode 
that  destroyed  the  C3  symmetry  axis,  and  predicted  a  Jahn-Teller  energy 
lowering  of  0.2  kcal  mole-1. 2i  A  more  recent  calculation,  including  cor¬ 
relation  effects  and  considering  the  coupling  of  all  three  e-type  vibra¬ 
tional  modes  for  methoxy,  predicts  different  effects:23 

(a)  Mode-coupling  increases  the  magnitude  of  the  Jahn-Teller  energy 
reduction. 

(b)  The  ground-state  electronic  energy  of  the  distorted  molecule 

is  0.75  kcal  mole-!  less  than  the  minimum  energy  for  the  Cx-u 
molecule.  v 


In  addition  to  optimizing  the  energy  relative  to  the  structural 
parameters  for  methoxy,  the  structural  parameters  for  methanol  were  also 
determined.  The  calculations  to  optimize  energy  with  respect  to  struc¬ 
ture  were  performed  using  the  fourth-order  MBPT  with  all  double-excita¬ 
tions,  but  the  core  orbitals  were  not  permitted  to  correlate  with  the 
valence  electrons.  A  variety  of  experimentally  determined  structural 
parameters24-^  and  the  parameters  predicted  by  these  calculations  appear 
in  Table  1.  The  agreement  between  the  predicted  parameters  and  the  re- 
empix ically-derived  parameters  of  Gerry,  et  al.,2^  is  excellent. 
Predicted  structural  parameters  for  the  methoxy  radical  are  presented  in 

~22  ” 

J  H.  A.  Jahn  and  E.  Teller ,  Proa.  Roy.  Soc. ,  Ser.  A ,  161 3  220  (1937). 

23 

G.  D.'  Bent,  G.  F.  Adams,  G.  D.  Purvis  and  R.  J.  Bartlett,  "The  Elec¬ 
tronic  Structure  of  Methoxy  Radical.  Jahn-Teller  and  Spin-Orbit  Energy 
Lowering" ,  manuscript  in  preparation. 

24 

J  K.  Kimuro  and  M.  Kubo ,  "Structure  of  Dimethyl  Ether  and  Methyl  Alcohol" , 
J.  Chem.  Phys .  30,  (1959). 


J  T.  Nishikawa,  no  title  available,  J.  Phys.  Soc.  Japan  11,  781  (1956). 

26 

P.  Venkateswarlu  and  W.  Gordy,  "Methyl  Alcohol.  I.  Microwave  Structure" 
J.  Chem.  Phys.  23,  1195  (1955).  3 


R.  M.  Lees  and  J.  G.  Baker,  " Torsion-Vibration— Rotation  Interactions  in 
Methanol.  I.",  J.  Chem.  Phys.  48  5299  (1968). 


M.  L.  C.  Gerry,  R.  M.  Lees  and  G.  Winnewisser ,  "The  Torsion-Rotation 
Microwave  Spectrum  of  ^CRz^OR  and  the  Structure  of  Methanol ",  J.  Mol 
Spec.  61,  231  (1967). 
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TABLE  1.  STRUCTURAL  PARAMETERS  FOR  METHANOL,  CH_OH 
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Table  2.  After  determining  the  optimal  structures  for  methanol  and 
methoxy,  MBPT  calculations  including  the  correlation  of  the  core  orbital 
electrons  were  done. 


The  results  of  the  electronic  structure  calculations  for  methanol 
and  the  methoxy  radical  provide  part  of  the  information  needed  to  esti¬ 
mate  the  dissociation  energy  of  methanol.  In  addition  to  the  electronic 
part  of  this  heat-of-reaction,  there  are  contributions  from  the  vibra¬ 
tional,  rotational  and  translational  degrees  of  freedom.  Thus, 


ahr 


AuVib 

4HR  + 


AH 


rot 


+  AH 


tr 

R 


At  0  K,  rotational  and  translational  degrees  of  freedom  do  not  contribute 
to  AHr,  and  AH^1  equals  the  difference  between  the  zero-point  energies 
of  the  reactant  and  the  products.  Since  the  hydrogen  atom  has  no  vibra¬ 
tional  degrees  of  freedom,  only  methanol  and  methoxy  contribute  to  AHj^^5. 
Tables  3  and  4  contain  summaries  of  the  molecular  information,  including 
electronic  energy  predictions  and  vibrational  frequencies,  for  methanol 
and  methoxy,  respectively.  All  the  vibrational  frequencies  for  methanol 
arc  known  from  experimental  data,  but  only  three  methoxy  frequencies 
have  been  observed."-  Thus,  six  frequencies  for  methoxy  must  be  estimated. 
Estimation  of  the  heat-of-reaction  for  reaction  (1), 


CH3OH  ->  CH-0  +  H 


(1) 


is  outlined  in  Table  5. 

A  fourth-order  MBPT  calculation,  that  includes  all  single-,  double-, 
and  quadruple-excitation  diagrams  that  contribute  at  fourth-order,  pre¬ 
dicts  a  value  for  the  electronic  part  of  the  heat-of-reaction.  This 
calculation  predicts  that  the  dissociation  is  endothermic,  requiring  an 
energy  of  107.3  kcal  mole"1  to  break  the  oxygen -hydrogen  bond.  To  com¬ 
pare  this  number  with  experiment,  however,  the  correction  for  the  change 
in  zero-point  vibrational  energy  that  occurs  during  the  dissociation  must 
be  included.  This  correction  reduces  the  energy  required  for  dissocia¬ 
tion  at  0  K  to  98 . 3  kcal  mole-l.  This  number  corresponds  to  AH§,  . 

Since  most  experimentally  tabulated  values  for  the  heat-of-reaction  are 
given  for  T  =  300  k,  the  effects  of  translational  and  rotational  degrees 
of  freedom  on  the  dissociation  energy  must  be  included.  As  the  details 


F.  C.  Engelking ,  G.  B.  Ellison  and  W.  C.  Lineberger ,  "Laser  photode- 
uaohnent  spectrometry  of  methoxide ,  deuteromethoxide  and  thiomethoxide . 
Electron  affinities  and  vibrational  structure  of  CEJ93  CD70  and  CH  S" 

J.  Chem.  Phys.  69  3  1826  (1978).  S3  3s 
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TABLE  2.  STRUCTURAL  PARAMETERS  FOR  THE  METHOXY 
RADICAL,  CH3Oa 
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rCH' 
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rCH" 
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H'CH" 

108°  36 

H"CH" 

112°  6' 

H'CO 

1110  30 

H"CO 

103°  54 

aBond  length  units t  A. 


TABLE  3.  ELECTRONIC  ENERGY  PREDICTIONS  AND 
VIBRATIONAL  FREQUENCIES  FOR  METHANOL,  CH3OH 


Electronic 


Vibration 


Energy  (h) 

SCF 

-115.07437 

D-MBPT (4) 

-115.49415 

SD-MBPT(4) 

-115.49713 

SDQ-MBPT (4) 

-115.49199 

Mode 

Frequency  (cm  ^ ) 

OH  stretch 

3681 

CH3  stretch 

3005 

CH3  stretch 

2965 

CH3  stretch 

2844 

Methyl  bend 
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Methyl  bend 

1477 

Methyl  bend 

1455 

COH 

1345 

CO  stretch 

1073 

Methyl  rock 

1165 

Methyl  rock 

1060 

torsion 

270 

19 


Electronic 


Vibration 


TABLE  4.  ELECTRONIC  ENERGY  PREDICTIONS  AND 
VIBRATIONAL  FREQUENCIES  FOR  METHOXY  RADICAL,  CHjO 


UHF 

D-MBPT (4) 
SD-MBPT (4) 
SDQ-MBPT (4) 
Mode 


Energy  (h) 
-114.45319 
-114.82484 
-114.82802 
-114.82321 
Frequency  (cm 


V 

C-H  stretch  (a^) 
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V 

H  umbella  (a^) 
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C-0  stretch  (a^) 

1020 

V 

C-H  stretch  (e) 

3010 

V 

H  bend  (e) 

1380 

V 

rock  (e) 
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TABLE  5.  DETAILS  OF  THE  ESTIMATION  OF  THE  DISSOCIATION 
ENERGY  FOR  METHANOL 


Electronic  Energy  (hartree) 

(Cel 

AE 

Zero-point  vibrational  energy 
(kcal  mole“l) 


,vib 


AE 


AEr(0  K) 

Translational  Energy  (300  K) 
(kcal  mole-1) 

.  ctr 
AE 

Rotational  Energy  (300  K) 
(kcal  mole-1) 


rot 


AE 


AEr  (300  K) 

AH°  (300  K) 
K 


CH3OH  -+ 
-115.49199 
. 17099h 
31.1 


ch3° 


H 


0.9 


0.9 


-114.82321  -.49778 

107.3  kcal  mole  * 

22.1  0. 


-1 


-9.0  kcal  mole 
98.3  kcal  mole-'*' 
0.9  0.9 

+0.9  kcal  mole-'* 
0.9  0. 


0. 


99.2  kcal  mole 


-1 


AE°  (300  K)  +  RT  =  99.8  kcal  mole 
K 


-l 


Experimental  values: 

100.9  kcal  mole  *  (Reference  29) 
104.0  kcal  mole  *  (Reference  30) 
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in  Table  5  show,  the  translational  energy  increases  during  the  course  of 
the  reaction,  while  the  rotational  energy  does  not  change.  Summing  the 
various  energy  contributions,  the  theoretical  analysis  predicts  that  the 
energy  of  reaction  is  99,2  kcal  mole-!.  Comparison  with  observed  enthal¬ 
pies  requires  addition  of  a  term  A(pV)  to  the  energy.  This  term  is  eval¬ 
uated  by  using  the  expression,  A(pV)  =  An(RT).  For  T  =  300  K,  the  A(pV) 
term  contributes  0.6  kcal  mole"-*-,  giving  an  estimated  dissociation  energy, 
AH^,3qo  =  99.8  kcal  mole-!. 

Two  experimental  values  for  the  dissociation  energy  of  methanol  have 
been  reported  recently.  Batt  and  McCulloch  derived  AH$,300  for  methanol 
using  thermochemical  parameters  derived  from  chemical  kinetic  data. 30 
They  derived  a  value  of  3.9  kcal  mole”!  for  heat-of-formation  of  the 
methoxy  radical,  leading  to  a  predicted  methanol  dissociation  energy  of 
104.0  kcal  mole-!.  Engelking,  et  al.,  measured  the  electron  affinity  of 
the  methoxy  radical,  and  estimated  that  the  heat-of-formation  for  CH3O 
was  0.7  kcal  mole-!. 29  This  heat-of-formation  required  that  the  methanol 
dissociation  energy  equal  100.9  kcal  mole-!.  The  results  of  the  MBPT  cal¬ 
culations  support  the  experimental  value  reported  by  Engleking,  et  al. 

If  the  heats-of-formation  for  methanol  and  hydrogen  are  assumed  to 
be  -48.0  kcal  mole-!  and  52.1  kcal  mole-!,  respectively,  the  theoretical 
calculations  predict  that  the  heat-of-formation  of  the  methoxy  radical  is 
0.3  kcal  mole-!.  A  summary  of  the  calculation  to  estimate  this  quantity 
appears  in  Table  6.  This  estimate  agrees  well  with  the  experimental 
value  of  Engelking,  et  al.29  This  prediction  can  be  checked,  however, 
by  calculating  the  dissociation  energy  for  the  methoxy  radical  and  com¬ 
bining  that  result  with  the  known  heats-of-formation  of  formaldehyde 
and  hydrogen.  This  calculation  is  presented  in  the  next  section. 

In  general,  good  quality,  correlated,  ab  initio  calculations  tend  to 
underestimate  heats-of-reaction.  Even  the  excellent  calculations  of 
Goddard  and  Schaefer  underestimated  the  formaldehyde  dissociation  energy 
by  about  7  kcal  mole-!./  This  characteristic  is  due,  in  part,  to  the 
finite  basis  sets  employed  in  molecular  calculations,  since  such  a  basis 
always  describes  the  products  better  than  it  describes  the  reactant. 

Such  a  condition  guarantees  that  the  electronic  energy  of  the  products  is 
relatively  lower  than  that  of  the  reactant.  Reliable  estimates  of  the 
magnitude  of  the  error  introduced  are  not  available,  but  the  fact  that 
theoretical  techniques  underestimate  dissociation  energies  is  general  and 
must  be  emphasized. 


30 

L.  Batt  and  R .  D .  McCulloch j  " Pyrolysis  of  Dimethyl  Peroxide ",  Int . 
J.  Chemical  Kinetics  8j  491  (1976). 

31 

G .  F.  Adorns^  nA  Priori  Estimation  of  Rate  Constants  for  Unimolecular 
Decomposition  Reactions” ^  ARBRL-TR-02143  (1979).  (AD#A069132) 
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TABLE  6. 

ESTIMATE  OF  METHOXY  RADICAL  HEAT-OF- FORMATION 

ch3oh 

CH30  +  H 

AHR,300 

99.8  kcal  mole  1 

AHf ,300 

(cH3°)  =  ah°j300  -  ah»j300  (H)  +  ah°j30Q  (ch3oh) 

AHf, 300 

(H)  .  '  52.1  kcal  mole  1 

AH£,300 

(CH3OH)  -48.0  kcal  mole"1 

"SSo 

(CFLO)  -0.3  kcal  mole-1 
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V.  DISSOCIATION  OF  THE  METHOXY  RADICAL: 

CH30  -*■  CH20  +  H 

The  dissociation  of  the  methoxy  radical  to  form  formaldehyde,  CH2O, 
and  hydrogen  is  a  reaction  that  appears  in  almost  all  chemical  models 
that  describe  the  oxidation  of  formaldehyde,  or  that  contain  significant 
contributions  by  formaldehyde  chemistry.  Experimental  studies  of  the 
dissociation  reaction  are  hampered  by  the  lack  of  molecular  information. 
A  theoretical  study  of  the  radical  offers  opportunities  to  estimate  the 
dissociation  energy  and  heat-of-formation  of  the  radical,  and  to  predict 
a  value  of  the  critical  energy  that  is  required  for  dissociation. 


Molecular  information  characterizing  the  methoxy  radical  was  dis¬ 
cussed  in  the  previous  section,  and  the  data  are  summarized  in  Table  4. 
fo  predict  the  dissociation  energy  for  the  methoxy  radical  requires  elec¬ 
tronic  energy  and  vibrational  frequency  values  for  formaldehyde.  Formal¬ 
dehyde,  CH2O,  is  a  closed-shell  molecule  with  spatial  symmetry  corre¬ 
sponding  to  the  point  group,  C2V.  Many  publications  describe  prior 
theoretical  studies  of  the  formaldehyde  molecule,  with  the  most  important 
of  these  summarized  in  the  recent  paper  by  Goddard  and  Schaefer.7 

Equilibrium  molecular  structure  parameters  for  formaldehyde,  sum¬ 
marized  in  Table  7,  were  predicted  by  optimizing  the  energy  as  a  function 
of  three  parameters;  the  CH  bond  distance,  the  CO  bond  distance  and  the 
HCH  bond  angle,  using  frozen-core  D-MBPT(4)  calculations.  Excellent 
agreement  obtains  between  the  predicted  and  the  experimental  values  for 
the  structural  parameters .32  Electronic  energy  predictions  and  vibra¬ 
tional  frequencies32  for  formaldehyde  appear  in  Table  8.  The  electronic 
energy  predicted  by  the  fourth-order  MBPT  calculation  that  includes  all 
single-,  double-  and  quadruple-excitation  diagrams  that  contribute  at 
that  order  is  lower  than  that  predicted  by  any  of  the  previous  theoretical 
studies  of  formaldehyde.  However,  because  of  the  non-variational  property 
of  the  1 inked-diagram  related  calculation,  this  energy  cannot  be  identi¬ 
fied  as  an  upper  bound  to  the  true  energy. 

Details  for  the  dissociation  energy  calculation  are  presented  in 
Table  9.  Theory  predicts  a  significantly  less  endoergic  dissociation 
energy  for  the  methoxy  radical  than  for  formaldehyde.  On  the  other  hand, 
the  predicted  dissociation  energy,  22.2  kcal  mole"1,  is  likely  less  than’ 
the  critical  energy  for  dissociation,  since  it  is  expected  that  the  dis¬ 
sociation  of  methoxy  will  have  an  energy  barrier  to  overcome.  In  this 
way  the  reaction  resembles  the  dissociation  reaction  of  the  formyl  rad¬ 
ical  . 2U 

Having  predicted  a  value  for  the  heat-of-reaction,  an  estimate  of 
the  heat-of-formation  for  the  methoxy  radical  can  be  obtained.  The  arith¬ 
metic  of  this  estimate  appear  in  Table  10.  With  the  heat-of-formation 


K.  Yamada,  T.  Nagakura^  K.  Kuchistu  and  Y.  Morino }  " Band  Contour 
Analysis  of  the  vj  and  x>2  Fundamentals  of  Formaldehyde J.  Mol. 
Spectrosoop.  38s  70  (1971). 
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TABLE  7. 

STRUCTURAL  PARAMETERS  FOR 

FORMALDEHYDE3 

Reference  32 

This  Work 

rch 

1.099 

1.102 

70 

n 

0 

1.203 

1.211 

HCH 

116°30 1 

116°11 ' 

aBond  length  units 3  A. 


TABLE  8.  ELECTRONIC  ENERGY  PREDICTIONS  AND  VIBRATIONAL 
FREQUENCIES  FOR  FORMALDEHYDE,  CH2O 


Electronic 


Vibration 


Energy  (h) 


SCF 

D-MBPT (4) 

SD-MBPT(4) 

SDQ-MBPT(4) 

Mode 

CH  stretch  (b2) 
CH  stretch  (a^) 
CO  stretch  (ai) 
CH2  bend  (ai) 
CH2  bend  (b2) 
CH2  rock  (b  ) 


-113.8974 
-114.28321 
-114.28868 
-114.28296 
Frequency  (cm 
2843 
2766 
1746 
1501 
1251 
1167 


) 
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TABLE  9.  DETAILS  OF  THE  ESTIMATION  OF  THE  DISSOCIATION 
ENERGY  FOR  METHOXY  RADICAL 


CH^O  ■+■ 

ch2o  + 

H 

Electronic  Energy 
(hartree) 

-114.82321 

-114.28296 

-.49778 

AEel 

•04238h  =  26.5  kcal 

mole  * 

Zero-point  vibrational 
energy  (kcal  mole-1) 

22.1 

16.2 

0. 

AEvib 

-5.9  kcal  mole"^ 

AEr(0  K) 

20.7  kcal  mole  * 

Translational  Energy 
(300  K0)  (kcal  mole-1) 

0.9 

0.9 

0.9 

4Etr 

0.9  kcal  mole  ^ 

Rotational  Energy 
(300  K) 

0.9 

0.9 

AErot 

0.0 

AE°  (300 

K) 

21.6  kcal  mole-^ 

AH° 

R,300 

E°  (300  k)  +  RT  = 

22.2  kcal  mole-^ 
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TABLE  10.  ESTIMATE  OF  THE  METHOXY  RADICAL  HEAT-OF -FORMAT ION 


CH30  + 

CH20  +  H 

AHR,300 

22.2  kcal  mole  ^ 

(CH30)  =  AHJ,  300^2°) 

+  AHf,300^  -  AHR,300 

AHf,300^CH20^ 

-26.0  kcal  mole  ^ 

AH£,300W 

52.1  kcal  mole  ^ 

AHf^300^CH30^ 

3.9  kcal  mole”'*' 
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of  formaldehyde  assumed  to  be  -26.0  kcal  mole-!,2’33  the  heat-of-forma- 
tion  result  is  3.9  kcal  mole-!,  a  value  that  differs  by  about  4  kcal 
mole-  from  the  value  obtained  using  the  methanol  dissociation  energy. 

A  major  part  of  the  difference  between  the  two  estimates  results  from 
the  systematic  underestimation  of  dissociation  energies  by  theoretical 
methods.  Note  that  the  error  decreases  the  heat-of-formation  estimate 
obtained  from  the  methoxy  radical  dissociation  energy.  More  accurate 
predictions  of  AHf  require  an  estimate  of  the  magnitude  of  the  error  in 
predicted  dissociation  energies. 


VI.  DISSOCIATION  OF  FORMALDEHYDE  TO  RADICAL  PRODUCTS: 

CH20  -»•  CHO  +  H 

Thermal  dissociation  of  formaldehyde  produces  radical  products  hydro¬ 
gen  atom  (2S)  and  formyl  radical,  HCO  (X2A').  Interestingly,  the  dis¬ 
sociation  to  molecular  products  H2(Ie+)  and  C0(Ie+)  is  almost  thermo¬ 
neutral,  while  the  observed  dissociation  energy  for  formaldehyde  is 
approximately  88  kcal  mole-!. 33  Both  sets  of  products  correlate  by 
symmetry  with  the  ground  electronic  state  of  formaldehyde,  S0(XlAi). 
Electronic  excitation  of  formaldehyde  to  the  first  singlet  excited  state, 
Si (a!A' ) ,  requires  less  energy  than  the  threshold  energy  for  dissociation. 
Photochemical  experiments  indicate  that  excitation  of  formaldehyde  to  the 
Si  state  leads  to  three  major  processes, 34 


H2C0(S1V')  H0  +  CO 


(!) 


-»■  H  +  HCO 


(II) 


-*■  H_C0 (S  ,V")  +  hv. 
2  o 


(III) 


Experimentalists  generally  agree  that  process  (I) ,  production  of  molecu¬ 
lar  products,  dominates  for  energies  near  the  S\  origin  (3.495  eV,  80.6 
kcal  mole-1).  The  importance  of  the  radical  process  increases  with  in¬ 
creasing  excitation  energy.  For  an  excitation  energy  of  91.3  kcal  mole-! 
the  sum  of  quantum  yields  for  processes  (I)  and  (II)  is  almost  one  with 
molecular  quantum  yield  0.32  ±  0.03  and  the  radical  quantum  yield 


33 

P.  Warneck3  "Photoionization  von  Methanol  und  Formaldehyd" 3  Z. 
Naturforsoh  A263  2047  (1971). 

34 

R.  D.  MoQuigg  and  J.  G.  Calvert 3  "The  Photodeaomposition  of  CE203  CD203 
CHDO  and  CH2O-CD2O  Mixtures  at  Xenon  Flash  Lamp"3  J.  Am.  Chem.  Soo. 

913  1590  (1969). 
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0.68  ±  0.03.35  Although  many  experimental  studies  of  the  photodissoci¬ 
ation  of  formaldehyde  have  been  reported,  there  is  little  agreement  as 
to  the  detailed  mechanism  that  would  explain  the  observed  chemistry. 

Goddard  and  Schaefer  reviewed  the  various  mechanistic  possibilities 
in  their  report  on  extensive  configuration  interaction  (Cl)  studies  of 
the  features  of  the  potential  energy  surfaces  relevant  to  the  photodis¬ 
sociation  of  formaldehyde . 7  These  include  the  dissociation  to  radical 
products  by  breaking  a  CH  bond,  the  dissociation  to  molecular  products 
via  an  intermediate  on  the  SQ  surface,  and  the  rearrangement  of  formal¬ 
dehyde  to  hydroxycarbene,  HCOH,  on  the  SQ  surface.  The  Cl  study  pro¬ 
vided  improved  values  for  energy  barriers  for  dynamical  processes  in 
formaldehyde  photochemistry,  including  the  prediction  that  the  energy 
required  for  photochemical  isomerization  to  hydroxycarbene  is  only  slightly 
greater  than  the  energy  needed  for  dissociation  to  either  radical  or 
molecular  products.  Mechanistic  studies  of  formaldehyde  photodissociation 
should,  therefore,  consider  the  presence  of  the  hydroxycarbene  molecule. 

The  Cl  calculations  also  predict  small  barriers  above  the  Si  band  origin 
to  any  photochemical  reaction  of  formaldehyde. 

The  Cl  calculation  for  the  dissociation  to  radical  products  predicts 
the  heat-of-reaction  to  be  about  79  kcal  mole-1,  a  value  below  the  experi¬ 
mental  threshold  of  about  86  kcal  mole-1.36  ci  calculations,  in  general, 
underestimate  dissociation  energies  and  overestimate  barriers  to  reaction, 
in  part  due  to  variations  in  the  MO  basis  as  the  molecular  structure  is 
changed.  Many-body  perturbation  theory  uses  the  same  diagrams  at  every 
point  on  the  potential  energy  surface,  a  fact  that  may  reduce  the  in¬ 
herent  error  in  a  theoretical  study.  The  dissociation  of  formaldehyde 
to  radical  products  provides  an  excellent  example  to  test  this  contention. 

The  results  of  MBPT  calculations  for  the  ground  electronic  state  of 
formaldehyde  were  presented  in  Table  8.  Electronic  structure  calculations 
for  the  formyl  radical  include  a  thorough  study  of  structures  surrounding 
the  equilibrium  conformation,  as  well  as  a  number  of  structures  in  the 
hypersurface  associated  with  dissociation  to  hydrogen  and  carbon  monoxide. 
The  predicted  structural  paramenters  for  the  equilibrium  conformation 
are  given  in  Table  11,  along  with  parameters  deduced  from  a  number  of 


35 

A .  Horowitz  and  J.  G.  Calvert,  " The  Quantum  Efficiency  of  the  Primary 
Processes  in  Formaldehyde  Photolysis  at  3130A  and  25%Crt,  Int.  J.  Chem. 
Kinet.  10,  713  (1978). 
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J.  H.  Clark,  C.  B.  Moore  and  N.  S .  Nogar,  "The  Photochemistry  of  For¬ 
maldehyde.  Absolute  Quantum  Yields,  Radical  Reaction  and  NO  Reactions ", 
J.  Chem.  Rhys.  68,  1264  (1978). 
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TABLE  11.  STRUCTURAL  PARAMETERS  FOR  THE  FORMYL  RADICAL,  HCQ: 
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experimental  studies .37-42  Because  of  the  more  precise  data  obtained 
by  Brown  and  Ramsey ,42  the  structure  predicted  in  that  publication  should 
be  more  accurate  than  the  other  empirical  structures.  The  theoretical 
bond  lengths  and  angle  predictions  agree  well  with  the  Brown  and  Ramsey 
parameters . 

A  summary  of  the  electronic  structure  calculations  for  the  predicted 
equilibrium  conformation  of  the  radical  appears  in  Table  12.  The  wave 
function,  not  a  pure  spin  state,  has  a  multiplicity  equal  to  2.01,  indi¬ 
cating  that  the  zero-order  wave  function  contains  only  small  contamina¬ 
tion  from  excited  states.  The  vibrational  frequencies  listed  in  Table 
12,  taken  from  the  data  of  Brown  and  Ramsey,  correspond  to  a  vibrational 
zero-point  energy  of  7.8  kcal  mole-1 * * * S. *  for  the  radical. 

The  heat-of-reaction  estimate  for  the  dissociation  of  formaldehyde 
to  radical  products,  outlined  in  Table  13,  parallels  the  previous  dis¬ 
sociation  energy  calculations.  An  electronic  energy  barrier  of  93.4  kcal 
mole-1,  predicted  by  the  MBPT  calculations,  is  significantly  higher  than 
the  Cl  barrier,  88.5  kcal  mole-1.?  Since  the  zero-point  energy  differ¬ 
ences  are  virtually  identical,  the  predicted  dissociation  energy  for  the 
threshold  reaction,  85  kcal  mole-1,  is  larger  than  the  value  predicted 
by  the  Cl  calculations,  80.2  kcal  mole-1.7  As  noted  by  Goddard  and 
Schaefer,  the  Cl  calculations  underestimate  the  threshold  by  about 
5  kcal  mole-1. 


37 

G.  Herzberg  and  D .  Ramsey ^  " The  7500  to  4500  A  Absorption  System  of  the 
Free  HCO  Radical",  Proa.  Roy .  Soo .  (London)  A233 J  34  (1955). 
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Acta  A233  737  (1969). 
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HCO  and  DCO  Free  Radicals  by  Microwave  Spectroscopy" 3  Astrophys.  J. 

178 3  L95  (1972). 
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TABLE  12. 

ELECTRONIC  ENERGY  PREDICTIONS  AND  VIBRATIONAL  FREQUENCIES 
FOR  FORMYL  RADICAL,  HCO 

Electronic 

Energy  (h) 

UHF 

-113.2769 

D-MBPT(4) 

-113.6343 

SD-MBPT(4) 

-113.64155 

SDQ-MBPT(4) 

-113.63634 

Vibration 

Mode 

Frequency (cm  j 

CH  Stretch 

2488 

CO  Stretch 

1861 

HCO  Bend 

1080 
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TABLE  13.  DETAILS  OF  THE  ESTIMATION  OF  THE  DISSOCIATION 
ENERGY  FOR  FORMALDEHYDE 


CH20 

CHO  + 

H 

Electronic  Energy 
(hartree) 

-114.28296 

-113.63634 

-.49778 

,rel 

AE 

. 14884h  = 

93.4  kcal 

Zero-point  Vibrational 
Energy  (kcal  mole-1) 

16.2 

7.8 

0. 

AOVib 

AE 

-8.4  kcal  mole 

-1 

AEr(0  K) 

85.0  kcal  mole 

-1 

Translational  Energy 
(kcal  mole"l) 

0.9 

0.9 

0.9 

AEtr 

0.9  kcal  mole 

-1 

Rotational  Energy 
(kcal  mole~l) 

0.9 

0.9 

0. 

Acr0t 

AE 

0. 

AEr  +  pV  = 

86.5 

kcal 

mole  1 

Experimental  Values,  AH° 

o 

o 

85 

-  87 

kcal  mole  ^ 

(36) 

86.9 

±  2 

kcal  mole  ^ 

(43) 

88.2 

±  1. 

6  kcal  mole  ^ 

(33) 

33 


The  MBPT  calculations  predict  a  dissociation  energy  slightly  below 
the  experimental  values  of  Walshand  Benson43  and  Warneck.33 


The  predicted  dissociation  energy  leads  to  an  estimate  of  the  heat- 
of-formation  of  the  formyl  radical  by  a  calculation  outlined  in  Table 
14.  The  predicted  value,  AHf^ootHCO)  =  8.4  kcal  mole-!,  is  slightly 
less  than  the  value  inferred  from  Warneck 's  recent  study,  9.9  kcal 
mole-!. 33  As  noted  previously,  underestimation  of  the  heat-of-reaction 
leads  to  underestimation  of  the  heat-of-formation  of  a  product  species, 
so  this  result  does  not  surprise. 


Before  discussing  alternate  dissociation  mechanisms  for  formaldehyde, 
the  dissociation  of  formyl  radical  will  be  investigated,  in  order  to  ob¬ 
tain  another  estimate  of  AHf>300  for  the  radical. 

VII.  DISSOCIATION  OF  THE  FORMYL  RADICAL:  CHO  — *•  CO  +  H 

A  recent  paper  by  Adams,  et  al.,  describes  a  complete  theoretical 
study  of  the  formyl  radical  using  linked-diagram  related  methods. 20  This 
section  reviews  the  results  of  that  publication. 

Formyl  radical,  HCO,  dissociates  into  ground  state  hydrogen  and  car¬ 
bon  monoxide  along  a  single  potential  energy  surface.  This  surface  is 
expected  to  have  a  substantial  energy  barrier  (>  5  kcal  mole-!)  above 
the  dissociation  energy  of  the  radical.  The  major  objectives  of  the 
theoretical  study  of  the  formyl  radical  potential  energy  hypersurface 
were  to  predict  the  height  of  this  energy  barrier  and  to  specify  the 
structure  of  the  activated  complex  at  the  barrier.  Table  15  lists  the 
predicted  structural  parameters  for  the  transition  state.  The  barrier 
height,  (the  difference  in  electronic  energy  between  the  equilibrium  con¬ 
figuration  and  the  saddle  point)  corresponds  to  a  critical  energy  of 
18.5  kcal  mole”  .  The  predicted  transition  state  properties  were  used 
in  an  a  priori  statistical  analysis  of  the  dissociation  kinetics  for  the 
radical.44 

To  estimate  the  dissociation  energy  of  the  formyl  radical,  the 
electronic  energy  of  carbon  monoxide  is  required.  Table  16  lists  the 
results  of  the  MBPT  calculations  as  well  as  the  value  assigned  to  the 
CO  stretching  frequency.  Optimization  of  the  CO  bond  length  resulted  in 
exact  agreement  between  the  experimental  and  theoretical  values,  2.1276b. 


R.  Walsh  and  S.  W.  Benson3  "Kinetics  and  Mechanisms  of  Gas-Phase 
Reactions  between  Iodine  and  Formaldehyde  and  the  Carbon-Hydrogen  Bond 
Strength  in  Formaldehyde" ,  J.  Am.  Chem.  Soc.  88 3  4570  (1966). 
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G.  F.  Adams 3  "Estimation  of  a  Rate  Coefficient  for  the  Decomposition  of 
the  Formyl  Radical" 3  manuscript  to  be  submitted. 
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TABLE  14.  ESTIMATE  OF  THE  FORMYL  RADICAL  HEAT -OF -FORMAT I ON 


AH 


£,300 


(HCO) 


ch2o 

AHR,300 


=  AH 


R,  300 


AH 


£,300 


AH£,300 


AH 


£,300 


-  AH 


£,300 


(H) 

(CH20) 

(HCO) 


+  CHO  +  H 

86.5  kcal  mole 

(")  *  4Hf,300  (CH2°> 

52.1  kcal  mole 

-26.0  kcal  mole 

8.4  kcal  mole 


TABLE  15.  MOLECULAR  PARAMETERS  OF  THE  FORMYL  RADICAL  AT  THE  SADDLEPOINT 


Electronic  Energy 
(Hartrees) 


UHF 

SD-MBPT(4) 
Barrier  Height 

Structural  Parameters 


C2A') 

-113.2488 

-113.6048 

23.1  kcal  mole 


R 


CH 

RCO 

HCO 


3.35  b 
2.15  b 
115° 


TABLE  16.  ELECTRONIC  ENERGY  PREDICTIONS  AND  VIBRATIONAL 
FREQUENCY  FOR  CARBON  MONOXIDE,  CO 


Electronic 


Energy  (Hartree) 


Vibrational 


SCF  -112.76646 
D-MBPT(4)  -113.11047 
SD-MBPT (4)  -113.11732 
SDQ-MBPT (4)  -113.11235 


Frequency 


,  -1, 

(cm  ) 


CO  Stretch  2170. 
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1  he  calculation  to  predict  the  formyl  radical  dissociation  energy  is  sum¬ 
marized  in  Table  17.  This  example  differs  from  the  previous  estimates 
of  dissociation  energy,  since  the  change  in  rotational  energy  for  the 
reaction  differs  from  zero.  The  predicted  heat-of-reaction  is  lower 
than  the  suggested  experimental  values.33, 45 

The  predicted  heat-of-formation  of  the  formyl  radical,  detailed  in 
Table  18,  exceeds  the  value  obtained  in  the  calculation  based  upon  the 
dissociation  energy  of  formaldehyde.  As  was  the  case  with  estimates  of 
AHf  for  methoxy,  the  systematic  underestimation  of  dissociation  energies 
leads  to  predicted  values  of  AH^  that  differ  significantly  for  the  radical 
as  product  or  reactant.  To  estimate  heats-of-formation  more  accurately 
requires  data  about  the  magnitude  of  the  error  in  the  predicted  heat-of- 
reaction.  One  approach,  used  initially  by  Bauschlicher  and  Shavitt,46 
systematically  improves  the  basis  set  until  the  quantity  calculated  does 
not  vary  with  the  basis  set  change.  Although  this  technique  may  work  for 
a  particular  problem,  employing  such  an  approach  for  all  theoretical 
studies  requires  computer  resources  not  currently  available. 

An  alternative,  pragmatic  approach  is  to  use  the  verified  differences 
between  experimental  and  theoretical  values  of  dissociation  energies  to 
estimate  the  degree  to  which  theory  underestimates  the  physical  quantity. 
For  the  reactions  discussed  in  this  report,  Warneck*s  experimental  dis¬ 
sociation  energy  for  formaldehyde  is  least  subject  to  revision.  The 
theoretical  prediction  underestimates  Warneck!s  measurements  by  1.7  kcal 
mole-!.  This  example,  combined  with  the  results  published  by  Redmon,  et 
al.,19  suggests  that  theory  underestimates  dissociation  energies  by  ap¬ 
proximately  1.5  kcal  mole-!,  and  that  addition  of  this  quantity  to  the 
original  theoretical  estimates  corrects  the  dissociation  energy  in  a 
reasonable  way.  The  original  and  modified  theoretical  estimates  are  com¬ 
pared  with  experimental  dissociation  energies  in  Table  19.  This  modifi¬ 
cation  improves  the  agreement  between  theory  and  experiment,  a  not  unex¬ 
pected  result  of  this  choice  of  correction  factor.  Modified  heat-of- 
formation  estimates  for  methoxy  and  formyl  are  outlined  in  Tables  20  and 
21,  respectively.  Use  of  the  Correction"  factor  reduces  the  difference 
between  the  two  AH?  values  for  methoxy  to  1.2  kcal  mole^1  (2.1  kcal  mole"1 
and  3.3  kcal  mole-1),  and  reduces  the  difference  between  the  formyl  radi¬ 
cal  estimates  to  1.4  kcal  mole"1  (9.9  kcal  mole-!  and  11.3  kcal  mole-1). 
These  heat-of-formation  estimates  fall  within  the  error  limits  of  the  ex¬ 
perimental  values,  and  thus  it  seems  justified  to  claim  that  this  method 
predicts  accurate  values  for  dissociation  energies  and  heats-of-formation. 

A  comment  about  the  intrinsic  accuracy  of  dissociation  energy  data 
clarifies  the  significance  of  the  theoretical  results.  Most  published 
dissociation  energies  have  error  limits  of  several  kilocalories  per  mole. 


M .  A.  Haney  and  J.  C.  Franklin a  "Excess  Energies  in  Mass  Spectra  of 
Some  Oxygen  Containing  Compounds",  Trans.  Far.  Soc.  65,  1794  (1969). 

C.  Bauschlicher  and  I.  Shavitt,  "Accurate  ab  initio  calculations  on 
Singlet-Triplet  Separation  in  Methylene" ,  J.  A.  Chem.  Soc.  100 ,  739 
(1978). 


36 


TABLE  17.  DETAILS  OF  THE  ESTIMATION  OF  THE  DISSOCIATION 


ENERGY  FOR  FORMYL 

RADICAL 

CHO 

CO  + 

H 

Electronic  Energy 
(hartree) 

-113.63634 

-113.11235 

-.49778 

,cel 

AE 

. 02621h  = 

16.4  kcal  mole  ^ 

Zero-point  Vibrational 

Energy  (kcal  mole-1) 

7.8 

3.1 

0. 

.  cvib 

AE 

-4.7  kcal  mole  * 

AEr(0  K) 

11.7  kcal  mole  ^ 

Translational  Energy 
(300  K)  (kcal  mole-1) 

0.9 

0.9 

0.9 

A  Ctr 

AE 

0.9  kcal 

mole 

Rotational  Energy 
(300  K)  (kcal  mole-1) 

0.9 

0.6 

0. 

,crot 

AE 

-0.3  kcal 

.  -1 
mole 

AE®  (300 

K) 

12.3  kcal 

mole  ■*■ 

t  .  n  a  _  n  _ _ ^  ^  _ 

-1 

AH°  goo  =  AE°  (300  K)  +  pV  =  12.9  kcal  mole-'*' 


Experimental  Values:  15.5  ±  1.5  Reference  33 

17  ±  2.  Reference  45 
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TABLE  18.  ESTIMATE  OF  THE  FORMYL  RADICAL  HEAT-OF-FORMATION 


CHO 


CO  +  H 


AHR,300 

“f,300  (CH0>  =  AHR,300  -  “f  ,300  tC0) 
“f  ,300 


12.9  kcal  mole 


-1 


AHf  ,300  ^ 


52.1  kcal  mole 


-1 


“boo  (“) 


-26.4  kcal  mole 


-1 


“f,300  'HC0> 


12.8  kcal  mole 


-1 
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TABLE  19.  MODIFIED  HEAT-OF-REACTION  PREDICTIONS  (kcal  mole 
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TABLE  20.  MODIFIED  HEAT-OF-FORMATION  ESTIMATES,  METHOXY  (kcal  mole"1) 
a.  CH3OH  +  CH30  +  H 


AHR,300 

101.3 

AHf ,300 

(CH3OH) 

O 

00 

1 

AH£,300 

(H) 

52.1 

(CH30) 

1.2 

b.  CH30  -►  CH20  +  H 


AHR, 300 

23.7 

AHf,300CH2°^ 

-26.0 

2.4 
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TABLE  21.  MODIFIED  HEAT-OF-FORMATION  ESTIMATES,  FORMYL  (kcal  mole  1 ) 
a.  CHo0  CHO  +  H 


AHR,300 

88.0 

AHf ,300 

(CH20) 

-26.0 

AHf ,300 

(H) 

52.1 

AHf’,300 

(CHO) 

9.9 

AHR, 300 

14.4 

AHf ,300 

(CO) 

-26.4 

11.3 
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For  molecules  such  as  formaldehyde  which  permit  sophisticated  experimental 
studies,  thermochemical  data  are  accurate  and  precise.  But  for  species 
such  as  the  formyl  radical,  often  the  suggested  empirical  values  lie  out¬ 
side  the  respective  suggested  error  limits.  The  error  limits  suggested 
by  Warneck33  for  the  dissociation  energy  of  formyl  barely  contain  the 
dissociation  energy  reported  by  Haney  and  Franklin. 45  The  MBPT  calcula¬ 
tions  predict  that  both  values  are  too  high,  and  that  this  dissociation 
energy  is  less  than  15  kcal  mole-!.  Since  both  theory  and  experiment 
have  minimum  error  limits  of  ±  3  kcal  mole-!,  the  difference  between  our 
theoretical  estimates  and  available  experimental  data  is  inconsequential. 

VIII.  DISSOCIATION  OF  FORMALDEHYDE  TO  MOLECULAR  PRODUCTS: 

ch2o  h2  +  CO 

The  formation  of  molecular  products  in  the  dissociation  of  formalde¬ 
hyde  occurs  via  one  of  several  competitive  mechanisms.  Two  possible  reactions 
channels  were  studied  by  Goddard  and  Schaefer,  leading  to  identification 
of  two  reaction  intermediates.?  One  intermediate  occurs  on  the  So 
potential  energy  curve  that  leads  to  direct  formation  of  the  molecular 
products,  H2  +  CO.  The  other  intermediate  occurs  on  the  S0  surface  be¬ 
tween  ground  state  formaldehyde  and  the  metastable  species,  hydroxycarbene . 

In  addition  to  determining  the  structure  of  each  intermediate,  Goddard 
and  Schaefer  predicted  the  vibrational  frequencies  for  the  two  evanescent 
species,  thus  facilitating  the  task  of  estimating  zero-point  energy  dif¬ 
ferences.  The  rearrangement  to  hydroxycarbene  provides  an  example  of  the 
importance  of  the  1,2-hydrogen  shift  in  chemistry. 

The  configuration  interaction  calculations  used  by  Goddard  and 
Schaeffer  underestimate  the  dissociation  energy  of  formaldehyde  (to  form 
radical  products)  by  about  5  kcal  mole"!,  an  error  considerably  larger 
than  that  attained  by  the  MBPT  calculation  reported  in  section  VI.  It 
is  important  to  investigate  the  predictions  of  MBPT  for  the  transition 
state  geometries  to  determine  whether  the  two  theoretical  approaches 
again  obtain  different  results. 

The  approach  to  these  MBPT  calculations  differs  from  that  reported 
in  the  previous  sections.  Rather  than  determining  the  optimal  structural 
parameters  for  the  transition  states,  this  study  uses  the  structures  re¬ 
ported  by  Goddard  and  Schaefer.  The  results  of  the  MBPT  calcuations  for 
the  molecular  products  and  rearrangement  transition  states,  as  well  as 
the  Goddard-Schaefer  predictions  for  the  real  vibrational  eigenvalues, 
are  presented  in  Tables  22  and  23,  respectively.  In  addition,  Table  24 
contains  the  results  of  a  "supermolecule"  calculation  for  the  molecular 
products,  H2  and  CO.  Tables  22,  23  and  24  contain  all  the  molecular 
data  needed  to  estimate  thermochemical  parameters  relevant  to  the  dis¬ 
sociation  of  formaldehyde  to  molecular  products. 
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TABLE  22.  ELECTRONIC  ENERGY  PREDICTIONS  AND  VIBRATIONAL 
FREQUENCIES  FOR  THE  MOLECULAR  PRODUCTS  TRANSITION  STATE 


Electronic 


Vibrational 


TABLE  23. 


Electronic 


Vibrational 


Energy  (hartree] 


SCF 

-113.73169 

D-MBPT(4) 

-114.13206 

SD-MBPT(4) 

-114.13991 

SDQ-MBPT (4) 

-114.13256 

Mode 

Frequency 

CH  Stretch 

2760 

CO  Stretch 

1654 

HCO 

1137 

HCO 

800 

HCH 

592 

ELECTRONIC  ENERGY  PREDICTIONS  AND  VIBRATIONAL  FREQUENCIES 
FOR  THE  HYDROXYCARBENE  TRANSITION  STATE 


SCF 

D-MBPT(4) 
SD-MBPT (4) 
SDQ-MBPT(4) 

Mode 

"out -of -plane" 

CH 

CH 

HCH 

HCH 


Energy  (hartree) 

-113.73666 

-114.13652 

-114.14547 

-114.13911 

Frequency  (cm  ^ ) 

3675 

2803 

2339 

1333 

1038 
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TABLE  24.  ELECTRONIC  ENERGY  PREDICTIONS  AND  VIBRATIONAL  FREQUENCIES 

FOR  MOLECULAR  PRODUCTS,  H2  +  CO 


Electronic 

Energy  (hartree) 

SCF 

-113.89766 

D-MBPT (4) 

-114.27622 

SD-MBPT (4) 

-114.28317 

SDQ-MBPT (4) 

-114.27797 

Vibrational 

Mode 

Frequency  (cm  j 

H,  Stretch 

4150 

CO  Stretch 

2150 
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Isomerization  of  formaldehyde  to  form  hydroxycarbene  exemplifies  a 
1,2-hydrogen  shift.  This  class  of  chemical  reactions  plays  an  important 
role  in  many  organic  reaction  mechnisms.  A  recent  review  by  Schaefer 
elucidates  the  main  features  of  the  1,2-hydrogen  shift,  while  providing 
a  survey  of  the  experimental  literature. ^  Goddard  and  Schaefer  located 
the  transition  state  for  the  isomerization  using  a  gradient-SCF  tech¬ 
nique.  7  This  calculation  also  provided  estimates  of  the  vibrational  force 
constants,  from  which  the  vibrational  frequencies  for  the  transition  state 
can  be  derived.  Interestingly,  the  zero-point  energy  of  the  transition 
state  exceeds  that  of  the  equilibrium  molecule  by  0.5  kcal  mole-!.  The 
estimate  of  the  barrier  height,  outlined  in  Table  25,  reflects  the  in¬ 
crease  of  the  critical  energy  relative  to  the  electronic  energy  barrier. 
The  SDQ-MBPT(4)  results  predict  an  electronic  energy  barrier  equal  to  90.3 
kcal  mole-!,  corresponding  to  a  critical  energy  equal  to  90.8  kcal  mole-!. 
This  number  exceeds  the  Goddard-Schaefer  prediction  for  the  barrier,  89.3 
kcal  mole-!,  and  the  predicted  heat-of-reaction  to  radical  products,  85 
kcal  mole-!.  However,  the  difference  between  the  two  linked-diagram 
method  calculations,  5.8  kcal  mole-!,  is  much  less  than  the  difference 
between  the  two  Cl  predictions,  9.3  kcal  mole-!. 

The  final  formaldehyde  dissociation  calculation,  outlined  in  Table 
26,  corresponds  to  the  direct  formation  of  molecular  products.  An 
electronic  energy  barrier,  94.4  kcal  mole-!,  is  reduced  by  the  zero- 
point  energy  change,  -6.6  kcal  mole-!,  to  yield  a  predicted  critical 
energy  for  this  dissociation  pathway  of  87.8  kcal  mole"!,  and,  like  the 
barrier  to  rearrangement,  exceeds  the  energy  required  to  dissociate  to 
radical  products. 

The  calculation  of  the  heat-of-reaction  for  the  dissociation  of 
formaldehyde  to  molecular  products  is  outlined  in  Table  27.  The  elec¬ 
tronic  energy  difference,  3.13  kcal  mole"!,  leads  to  a  predicted  heat-of- 
reaction,  -2.1  kcal  mole-!,  that  underestimates  the  suggested  empirical 
value  by  1.8  kcal  mole-!. 

The  results  reported  here  agree  with  the  Goddard-Schaefer  electronic 
structure  calculations  and  support  their  proposed  mechanisms  for  the 
photochemical  dissociation  of  formaldehyde  to  form  molecular  products. 

In  fact,  the  MBPT  results  predict  that  the  energy  requirements  for  the 
production  of  molecular  products  via  either  of  the  proposed  mechanisms 
approach  the  energy  requirements  for  thermal  dissociation  to  radical 
products.  To  determine  whether  these  results  are  an  artifact  of  the  MBPT 
method,  double-excitation  coupled-cluster  calculations  paralleling  the 
many-body  calculations  were  performed.  Table  28  contains  the  result 


47 

H .  F,  Schaefer,  III,  ", The  1,2-Hydrogen  Shift:  A  Common  Vehicle  for 
the  disappearance  of  Evanescent  Molecular  Species n,  Accts .  Chem. 
Research,  in  press . 
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TABLE  25.  BARRIER  HEIGHT  FOR  FORMALDEHYDE  TO  HYDROCARBENE  ISOMERIZATION 


h2co  -> 

Electronic  Energy  -114.28296 

(hartree) 


(HCOH)* 

-114.13911 


AE 


el 


. 14385h  =  90.3  kcal  mole 


Zero-point  Vibrational 
Energy*  (kcal  mole-1) 

AE  0.5  kcal  mole"'*' 

z,p 

AEo^  90.8  kcal  mole 


Zero-point  vibrational  energies  calculated  by  Goddard  and  Schaefer(7). 
Although  their  predicted  zero-point  energy  for  formaldehyde  is  greater 
than  the  experimental  value ,  the  change  in  zero-point  energy  between 
species  seems  to  be  insensitive  to  this  error. 


TABLE  26. 

BARRIER  HEIGHT  FOR  DISSOCIATION 

OF  FORMALDEHYDE  TO  MOLECULAR 

PRODUCTS 

ch2o 

(ch2o)* 

Electronic  Energy 

-114.28296 

-114.13256 

(hartree) 

AEel 

.1504h  =  94.4  kcal 

mole” 

Zero-point  Vibrational 
Energy*  (kcal  mole-!) 

18.3 

11.7 

AE. 

zp 

-6.6  kcal  mole  * 

'  AEo* 

87.8  kcal  mole"*' 

* 

Zero-point  energies  predicted  by  Goddard  and  Schaefer  (7). 
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TABLE  27.  HEAT- OF- REACT I ON  FOR  THE  DISSOCIATION 
OF  FORMALDEHYDE  TO  MOLECULAR  PRODUCTS 


Electronic  Energy 
(hartree) 


AE 


el 


Zero-point  Vibrational 
Energy  (kcal  mole-1) 


AE 

ae; 


Zp 


Translational  Energy 
(kcal  mole"l) 


AE 


tr 


Rotational  Energy 
(kcal  mole-1) 


AE 


rot 


ch20 


-114.28296 


H„ 


.00499h  =  3.13  kcal  mole 
16.2 


CO 

-114.27797 

9.1 


-7.1  kcal  mole 
-3.9  kcal  mole 
0.9  0.9 

0.9  kcal  mole 
0.9  0.6 

0.3  kcal  mole 
1 


-1 

-1 


0.9 


0.6 


AE°  (300  K)  -  2.7  kcal  mole 
AH°  300  =  E°  (300  K)  +  pV  =  -2.1  kcal  mole' 


TABLE  28.  RESULTS  OF  CCD  CALCULATIONS  FOR 
FORMALDEHYDE  CHEMISTRY 


Molecule 


CCD  Energy  (hartree) 


CH20,  equilibrium 

Transition  state,  rearrangement 
Transition  state,  molecular  products 

CHO,  equilibrium 


-114.27718 

-114.12930 

-114.12407 

-113.6288 
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of  the  CCD  calculations  of  the  electronic  energy  for  several  of  the 
molecular  and  intermediate  structures  relevant  to  formaldehyde  dissoci¬ 
ation,  while  Table  29  contains  a  comparison  of  the  Cl,  MBPT  and  CCD 
predictions  of  thermochemical  quantities.  The  CCD  calculations  predict 
the  highest  values  for  the  thermochemical  quantities,  with  a  notably 
large  barrier  predicted  for  the  rearrangement  reaction.  In  sum,  the  cal¬ 
culations  suggest  that  the  rearrangement  barrier  exceeds  the  barrier  for 
dissociation  to  molecular  products. 

In  their  discussion  of  the  Cl  results,  Goddard  and  Schaefer  comment 
that  configuration  interaction  calculations  generally  underestimate  heats- 
of-reaction  and  overestimate  barrier  heights.  The  expected  error  in  the 
dissociation  energy  calculations  equals  about  5  kcal  mole-1,  while  the 
expected  error  in  the  barrier  height  calculation  is  approximately  2  kcal 
mole-1.  Linked-diagram  calculations  have,  to  date,  underestimated  dis¬ 
sociation  energies  by  1  to  2  kcal  mole-!,  and  overestimated  rearrange¬ 
ment  barriers  by  similar  amounts. 19  Assuming  these  error  estimates,  the 
relative  energy  differences  indicated  in  Table  29  reduce  to  almost  zero. 
Therefore,  using  energy  as  the  criterion  for  reaction,  all  three  of  the 
proposed  reaction  channels  should  be  considered  when  seeking  a  mechanism 
to  describe  the  photo  dissociation  of  formaldehyde.  The  inferences  with 
respect  to  the  thermal  dissociation  of  formaldehyde  are  considered  in  the 
next  section. 


IX.  THE  UNIMOLECULAR  THERMAL  DISSOCIATION  OF  FORMALDEHYDE 

Combustion  chemistry  models  customarily  assume  that  the  thermal  dis¬ 
sociation  of  formaldehyde  produces  radical  products, 


CH20  CHO  +  H. 


Experimental  studies  by  Clark,  Moore  and  Nogar  demonstrated  that  the 
dissociation  to  radical  products  exhibits  an  abrupt  energy  threshold  be¬ 
tween  85  and  87  kcal  mole“1.36  Such  a  value  falls  within  the  error  bars 
of  the  thermodynamic  estimates  of  Walsh  and  Benson,  86.9  ±  2  kcal  mole-1, 
and  Warneck,  88.2  ±1.6  kcal  mole“1.33  Recent  work  by  Reilly,  Clark, 

Moore  and  Pimentel  confirms  the  radical  threshold  at  86  ±  1  kcal  mole“1.48 
The  MBPT  and  CCD  calculations,  corrected  for  T  =  300  K,  predict  the  rad¬ 
ical  threshold  at  86.6  and  87.8  kcal  mole-1,  respectively.  As  noted  in 
the  previous  section,  the  relative  energy  calculations  suggest  that  the 
dissociation  to  molecular  products  should  compete  with  the  radical  pro¬ 
duct  channel. 
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J •  P.  Reilly,  J.  H .  Clark,  C .  23.  Moore  and  G.  C .  Pimentel,  to  be  pub¬ 
lished.  Results  quoted  by  Goddard  and  Schaefer.? 
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TABLE  29.  COMPARISON  OF  Cl,  MBPT(4)  AND  CCD  CALCULATIONS 


(UNITS: 

kcal  mole-1) 

Reaction 

Cl 

SDQ-MBPT (4) 

CCD 

CH20  +  CHO  +  H 

^80. 

85.0 

86.2 

CH20  -»■  (HCOH)^ 

88.8 

90.8 

93.3 

ch2o  +  (CI^O)* 

87. 

87.8 

89.5 
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Several  other  results  of  the  photochemical  experiments  relate  to 
this  discussion.  These  experiments  involve  excitation  into  the  first 
excited  singlet  state,  Si,  and  subsequent  dissociation  to  either  radical 
or  molecular  products.  For  energies  below  the  radical  threshold,  formal¬ 
dehyde  (Si)  dissociates  to  molecular  products  with  quantum  yield  0.9. 

For  energies  equal  to  and  greater  than  the  radical  threshold,  the  total 
quantum  yield  is  near  1.0,  with  radicals  accounting  for  about  forty  per 
cent  of  the  dissociation.  The  abrupt  threshold  for  radical  products 
demonstrates  that  no  barrier  occurs  in  this  reaction  channel,  and  this 
implies  that  the  dissociation  does  not  occur  on  the  Si  hypersurface. 

No  mechanism  to  explain  the  photodissociation  of  formaldehyde  has 
yet  received  the  acceptance  of  the  chemical  community.  Goddard  and 
Schaefer f s  theoretical  study  suggests  that  formaldehyde  (Si)  undergoes 
internal  conversion  to  formaldehyde  (S0) ,  followed  by  dissociation  that 
is  enhanced  by  quantum  mechnical  tunneling  to  molecular  products.  These 
proposed  reactions  occur  on  the  SQ  hypersurface,  so  there  is  no  priori 
reason  why  molecular  products  may  not  occur  during  thermal  dissociation, 
since  the  transition  state  for  this  channel  is  energetically  competive 
with  the  dissociation  to  radical  products.  However,  dynamical  consider¬ 
ations  strongly  favor  the  formation  of  radical  products.  The  "loose” 
transition  state  representative  of  that  dissociaton  reaction  guarantees 
that  the  entropy  of  activation  for  this  channel  will  exceed  that  char¬ 
acteristic  of  a  "tight"  transition  state.  The  two  transition  states 
found  by  Goddard  and  Schaefer  are  "tight"  transition  states.  Consequently, 
the  dissociation  to  radical  products  should  occur  more  rapidly,  in  a 
chemical  kinetics  sense,  than  the  dissociation  to  molecular  products, 
provided  that  the  non-fixed  internal  energy  distribution  is  statistical . 

In  photodissociaton,  the  constancy  of  the  molecular/radical  branching 
ratio  above  the  radical  threshold  implies  that  the  reacting  species  does 
not  have  a  statistical  distribution  of  non-fixed  internal  energy, 

The  connection  between  the  confusion  associated  with  the  mechanism 
of  photodissociation  of  formaldehyde  and  the  simplistic  thermal  dis¬ 
sociation  mechansim  defies  easy  definition.  That  some  anomaly  in  thermal 
dissociation  may  occur  under  extreme  conditions  of  temperature  and  pres¬ 
sure  seems  likely,  but  it  may  be  that  under  extreme  conditions  the  dissoci¬ 
ation  reaction  is  ignorable  relative  to  competing  bimolecular  reactions. 


X.  SUMMARY 

The  calculations  described  in  this  report  have  important  implica¬ 
tions  for  research  in  flame  zone  chemistry.  For  the  first  time  an  a 
priori  model  predicts  heats-of-reaction  as  accurately  as  the  best  experi¬ 
ments.  In  addition,  heats-of-formation  can  be  estimated  with  moderate 
accuracy.  This  accuracy  is  improved  if  a  pragmatic  correction  to  the 
heat-of-reaction  estimate  is  used. 
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The  results  of  calculations  of  dissociation  energies  and  barrier 
heights  for  a  series  of  formaldehyde  reactions  support  the  photodissoci¬ 
ation  mechanism  proposed  by  Goddard  and  Schaefer. 7 

The  linked-diagram  based  calculations  described  in  this  report 
appear  to  predict  thermochemical  quantities  more  accurately  than  do 
truncated  configuration  interaction  calculations. 
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APPENDIX  A.  THE  CORRELATION  PROBLEM  AND  THE  METHOD 
OF  CONFIGURATION  INTERACTION 

The  standard  approach  to  electronic  structure  theory  for  many  re¬ 
searchers  is  the  Hartree-Fock  self-consistent  field  method,  a  technique 
based  upon  a  variational  principal.  For  a  normalized  approximate  wave- 
function  ip,  the  expectation  value  of  the  electronic  Hamiltonian  operator, 
corresponding  to  the  electronic  energy*  has  the  form, 


E  =  <>lHh>  (A- 1 ) 

The  variational  principle  guarantees  that  this  expectation  value,  E,  ex¬ 
ceeds  the  true  value  of  the  electronic  energy;  that  is,  the  expectation 
value  is  an  upper  bound. 

The  Hartree-Fock  wave-function  may  be  written  as  an  antisymmetrized 
product 


fHF  =  A(-m)  <i)i(1)(i)2^2^  •  • 


(A-2) 


where  A  is  the  antisymmetrizer  for  m-electrons  and  the  <j>ifs  are  one-elec¬ 
tron  spin  orbitals,  products  of  a  spatial  orbital  x  and  a  one-electron 
spin  function.  This  function  has  a  convenient  expression  as  a  Slater- 
determinant 


HF 


^(1)  *2(l)...+m(l) 

.  •  * 

•  •  .  . 

qO)  <j>2  (m) .  .  .<j>m(m) 


(A3) 


The  Hartree-Fock  wave-function  corresponds  to  the  determinant  for  which 
the  orbitals  in  equation  (2)  have  been  varied  to  give  the  lowest  possible 
electronic  energy . 

Minimizing  the  energy  resulting  from  a  single-determinant  wave- 
function  leads  to  a  complicated  set  of  integro-differential  equations, 
the  Hartree-Fock  equations.  The  conceptual  importance  of  the  Hartree- 
Fock  approach  is  that  the  HF  wave  function  is  the  best  wave  function 
that  can  be  constructed  by  assigning  each  electron  to  a  separate  orbital 
or  function  that  depends  only  upon  the  coordinates  of  that  electron. 

In  practice,  the  one  electron  functions,  are  expanded  in  terms 

of  a  finite  set  of  analytical  basis  functions. 
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CA-4) 


Because  of  the  finite  basis,  only  approximate  solutions  of  the  HF  equa¬ 
tions  occur.  The  lowest-energy  single-determinant  wave  function  avail¬ 
able  in  the  finite  basis,  the  so-called  self-consistent  field  (SCF)  wave- 
function,  represents  the  level  of  theoretical  sophistication  most  commonly 
encountered  in  the  chemical  literature.  SCF  wave-functions  suffice  to 
describe  the  structure  and  properties  of  molecules  near  equilibrium  in 
the  ground  electronic  state.  This  level  of  calculation  does  not  suffice 
to  describe  potential  energy  curves  and  surfaces,  transition  states,  or 
excited  electronic  states.  The  difference  between  the  best  Hartree-Fock 
approximation  of  the  electronic  energy  and  the  exact  energy  calculated 
using  the  non-relativistic  Schrodinger  equation  is  the  correlation  energy. 
The  name  reflects  the  physical  concept  that  the  Hartree-Fock  potential 
contains  the  average,  rather  than  the  instantaneous,  inter-electron 
potential,  and  thus  neglects  the  correlation  between  the  motions  of 
electrons.* 


To  circumvent  the  error  introduced  by  the  HF  approximation,  excited 
electronic  configurations  were  added  to  the  Hartree-Fock  wave  function, 
leading  to  the  method  of  configuration  interaction  (Cl).  A  Cl  wave 
function  can  be  represented  as  a  sum  of  Slater  determinants,  and  the 
energy  is  determined  by  variation  of  the  expansion  coefficients  until  a 
minimum  energy  is  obtained.  Were  the  basis  set  of  one  electron  functions 
(MO* s)  complete,  the  result  of  a  complete  Cl  calculation  would  approach 
the  exact  numerical  solution  of  the  N-particle  Schrodinger  equation. 

The  Cl  wave  function,  with  $0  the  HF  solution,  has  the  form. 


¥  =  $o  +  .  $a  + 

1  ,a  1  i 


„ab  .ab  , 

.  .  C .  .  $ .  .  + 
i>J  iJ  iJ 


(A-5) 


where,  for  instance,  $>.  represents  the  configuration  obtained  by  ex¬ 
citing  an  electron  from  orbital  i  of  $0  to  an  orbital  a  in  the  virtual 
orbital  space.  Equation  (A-5)  has  the  more  compact  form. 


4s  pointed  out  by  Shavitt  (1)3  this  definition  of  correlation  energy 
neglects  non-dynamical  effects  caused  by  near  degeneracies  and  re¬ 
arrangements  of  electrons  within  partially  filled  shells.  ■  This  leads , 
in  the  case  of  potential  energy  curves  that  are  always  poorly  described 
by  RHF  wave  functions 3  to  the  counterintuitive  idea  that  the  correlation 
energy  increases  as  the  electons  move  apart  with  the  separating  atoms. 

A1 

I.  Shavitt ,  "The  Method  of  Configuration  Interaction" ,  in  Modern 
Theoretical  Chemistry^  Vol_.  III.  Methods  of  Electronfa  Structure  Theorjjj 
edited  by  H.  F.  Schaefer 3  III3  Plenum3  Few  York3  1977. 


58 


(A-6) 


m 

Z 

s=l 


C  cj) 


where  the  $s  represent  pre-determined  expansion  functions  and  the  co¬ 
efficients  Cs  represent  the  variational  parameters  used  to  minimize 
E[y] .  This  form  of  the  wave-function  leads  to  the  generalized  matrix 
eigenvalue  equation, 


H  £  =  E  S_  c, 


(A-7) 


with  the  matrix  elements,  Hst  and  Sst,  defined  by 


Hst  5  <*slHl*t>  (A-8) 

Sst  =<4sl»t>.  CA-9) 


and  the  column  vector  c_  has  the  coefficients  cs  as  components. 

The  principal  computational  step  of  a  Cl  calculation  is  the  evalu¬ 
ation  of  the  matrix  elements,  Hst-  Choosing  the  molecular  orbital  basis 
set  to  be  orthonormal  facilitates  this  evaluation,  for  then  the  matrix 
elements  equations  have  the  form 


Hst  =<4>l 
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E 

i  J 


st, 

a. .h.  . 
13  13 
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ijkl 


,  st 

bijki 


’ijkl  ’ 


(A-10) 


with  the  orbital  integrals  defined  by 


(A— 11) 


where  g ^ 


r 


-1 

12* 


(A-12) 


The  orbital  integrals  are  obtained  from  analogous  basis-set  integrals 


K 

pq 


<Xp(rp 


hilxq 


(A- 13) 
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8pqrs  ^p^rl^r^r2^  ^q^rl^s  ^r2^’ 
by  means  of  the  transformation 

h.  .  =  Z  U*.  h  U  . 

ij  p>q  p1  pq  qj 


8ijkl  p,q^r,s  lJpi  Uqj  Urk  Usl  8pqrs, 


where 
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U  .. 
P  Pi 


(A-15) 


(A-16) 


(A-17) 


The  transformation  of  the  two-electron  integrals,  equation  (A-16),  requires 
extensive  computational  resources,  thus,  demanding  efficient  implementation 
of  this  step. 

The  computational  steps  required  for  a  Cl  calculation  can  be  outlined 
as  follows: 

1 .  Choose  an  analytical  basis  set  (xp)  and  compute  the  basis  set 
integrals,  equations  (A-13)  and  (A-14). 

2.  Choose  the  orbitals  (41),  by  an  SCF  or  similar  calculation,  and 
transform  the  basis-set  integrals  to  molecular  orbital  integrals. 

3.  Choose  and  construct  a  set  of  symmetry  adapted  configuration 
functions,  $,  appropriate  to  the  state  under  consideration,  and  compute 
the  Hamiltonian  matrix,  H,  in  terms  of  these  CF's. 

4.  Compute  the  several  lowest  eigenvalues  of  the  matrix  H. 

A  full  Cl  calculation,  including  all  possible  N-tuple  excitations 
for  the  N-electron  problem,  provides  the  exact  solution  to  the  problem. 

Full  Cl  calculations  satisfy  all  of  Pople's  desiderata. 2  Unfortunately, 
full  Cl  calculations  for  even  relatively  small  molecules  require  im¬ 
practical  amounts  of  computer  time  and  hardware  resources.  On  the  other 
hand,  a  truncated  Cl  expansion  suffers  from  three  interrelated  problems: 


A2 

J.  A.  Pople,  R.  irishman,  H.  B.  Schlegel  and  J.  S.  Binkley ,  "Electron 
Correlation  Theories  and  Their  Application  to  the  Study  of  Simple  Poten¬ 
tial  Surfaces ",  Int.  J.  Quantum  Chem. ,  14,  545  (1978). 
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1.  The  convergence  of  the  method  for  a  given  level  of  excitation 
is  slow. 

2.  The  number  of  configurations  increases  nonlinearly  with  the 
level  of  excitation  included. 

3.  Truncated  Cl  calculations  are  not  size-consistent. 

In  most  applications  of  Cl,  only  double-excitations  relative  to  a  single 
determinant  are  included  in  the  wave-function  expansion.  The  magnitude 
of  the  size-consistency  error  for  such  a  wave-function  may  be  reduced 
by  including  quadruple  excitations  in  the  wave -function,  or  by  replacing 
the  single  determinant  by  a  multi-determinental  reference  function, 
thereby  including  the  most  significant  quadruple  excitation  contributions. 

For  many  problems,  the  method  of  configuration  interactions  remains 
the  most  reliable  and  best  understood  theoretical  technique  available  for 
the  study  of  molecular  electronic  structure.  However,  the  difficulty  and 
expense  of  constructing  and  diagonalizing  the  Hamiltonian  matrix  have 
encouraged  alternative  methods  to  estimate  the  correlation  effect.  Each 
of  the  following  appendices  outlines  a  successful  alternative  to  Cl, 
while  future  reports  will  describe  the  development  of  a  powerful  new 
theoretical  method,  the  graphical  unitary  group  approach. 4»5 ,6 
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APPENDIX  B.  THE  COUPLED  CLUSTER  APPROACH 


The  coupled  cluster  approach  develops  an  alternative  to  the  Cl 
method  by  considering  the  origin  of  particular  levels  of  excitation  in 
CI.1>2  The  basic  element  in  CCA  is  a  cluster  expansion  where  one-body, 
two-body,...,  n-body  clusters  are  the  fundamental  entities.  A  coupled 
cluster  wave-function  is  represented  by. 


cc 


(B-l ) 


where  T  is  separated  into  the  various  cluster  terms. 


T=T1*T2*T3*  ... 


CB-2) 


In  the  occupation  number  representation. 


=.E  A+  A. 

i ,a  i  a  i 

(B-3) 

ah  .  +  + 

=  2  t.  .  A  A.  A.  A. 

a  1  b  J 

(B-4) 

a>b 


where  the  t's  are  anti symmetrized  components  to  be  determined.  Expanding 
the  exponential  operator  gives  a  new  expression  for  the  coupled  cluster 
wavefunction. 


J.  Cizek,  "On  the  Correlation  Problem  in  Atomic  and  Molecular  Systems. 
Calculation  of  Wave  Functions",  J.  Chem.  Phys  45,  4256  (1966). 

B2 

J.  Paldus  and  J.  Cizek,  "Relation  of  Coupled  Pair  Theory ,  Cl,  and  Some 
Other  Many -Body  Approaches",  Energy,  Structure  and  Reactiviti^,  Smith 
and  McCrae,  Eds.  (Wiley,  New  York,  1973). 
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*cc  =  (1  +  T1  +  T2  +  +  1/2  Tl  +  X/2  T2  +  +  T1  T2 


+ 


">  l»o> 


(B-5) 


Similary,  a  Cl  wavefunction  has  the  form. 


V 


Cl 


(1  +  C1  *  C2  . 


(E-6) 


where  the  C^s  refer  to  excitation-levels  relative  to  $  .  By  relating 
terms  in  equations  (B-5)  and  (B-6) ,  it  follows  that  0 


C2  -  T2  +  1/2  7* 

c3  =  t3  +  TxT2  +(1/3!)  Tj3 

C4  =  t4  +  1/2  T22  +  TjT^  +  fl/4!)^4  +  (1/21)7^7^ 


(B-7) 


The  quadruple-excitation  term  in  Cl  consists  of  several  different  parts: 
a  four -body  term,  a  product  of  two-body  terms  which  is  said  to  be  "dis¬ 
connected",  and  additional  disconnected  products  that  involve  Tl.  For 
3>o  represented  by  SCF-orbitals,  Brilloun's  theorem  justifies  ignoring 
all  terms  that  include  Tj  up  to  fourth-order.  Furthermore,  for  closed- 
shell  problems,  the  most  important  contributors  to  the  correlation 
energy  are  the  two-body  cluster  terms,  T2.  This  leads  to  the  basic  pre¬ 
mise  of  coupled-pair  theory,  that  T2  approximates  T. 2  in  particular, 
coupled-pair  many-electron  theory  (CPMET) ,  developed  by  Cizek,  uses  the 
wave-function, 

T 

^CPMET  =  e  l$o^=  t1  +  T2  +  1 /2  T2  +  (B_8) 


Pople  calls  this  approximation  coupled-cluster  double-excitation,  or 
CCD.  The  CCD  (CPMET)  wavefunction  has  the  explicit  form, 


B3 

V.  4.  Pople,  R.  Krishman,  H.  B.  Sohlegel  and  J.  S.  Binkley ,  "Electron 
Correlation  Theories  and  Their  Application  to  the  Study  of  Simple 
Reaction  Potential  Surfaces ",  Int.  J.  Quantum  Chem .  14,  545  (1978). 
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(B-9) 


CCD 


$ 

o 


+  1/4 


E  Aab  ab 
.  .  .  A.  .t.  . 
13  ab  13  13 


$ 

o 


+  1/32 


E  E  ab.cd  ab  cd 
ijab  klcd  ij  kl  ij  kl 


$ 

o 


The  operator  product  t^jt]^  is  zero  if  any  coincidences  occur  among  the 
eight  indices,  and  represents  a  quadruple  excitation  otherwise.  In  the 
latter  case,  T  becomes. 


CCD 


$  +  1/4 

o 


E 

ijab 


Aal?  +  1/32 

ij  ij 


E  E  .ab.cd  abed 
ijab  klcd  ij' kl  ijkl 


(B-10) 


Explicit  equations  for  the  coefficients,  A,  are  obtained  by  requiring 
that  the  projection  of  the  function 


<<H  -  ^bccD 


> 


C  B  —  1 1 ) 


ab 

on  $  and  V . .  equals  zero. 
0  ij  n 


Thus 


<*o  lH  -  EIW“  0 

<* i j  lH  -  El  W 


(B  - 1 2  ) 
(B —13) 


These  equations  suffice  to  determine  the  coefficients,  A,  and  the  energy, 
E.  Details  of  the  reduction  of  these  formal  equations  to  expressions 
including  integrals  involving  one-  or  two-electrons  are  given  by  Pople, 
et  al,3  and  Bartlett  and  Purvis. 4  Substitution  of  'I'cCD  (B-10),  in 
equation  (B-12)  yields  an  expression  for  the  energy  in  terms  of  two- 
electron  integrals  and  the  coefficients,  A.  Substitution  of  the  energy 
expression  in  C B- 13)  permits  elimination  of  energy  as  a  variable,  and 
leaves  an  equation  relating  the  coefficients,  A,  with  the  one-  and  two- 
electron  molecular  integrals.  That  equation,  however,  requires  an 
iterative  solution. 


B4 

R.  J.  Bartlett  and  G.  D.  Purvis 3  "Many-Body  Perturbation  Theory 3 
Coupled-Pair  Many-Eleotron  Theory 3  and  the  Importance  of  Quadruple 
Excitations  for  the  Correlation  Problem" 3  Int.  J.  Quantum  Chem.  14  3 
561  (1978). 
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As  implemented  by  Pople,  et  al,^  and  Bartlett  and  Purvis,**  the  CCD 
met  od  includes  excitations  relative  to  a  single-determinant  expression 
tor  $0.  The  technique  satisfies  the  size-consistency  and  invariance 
desiderata,  but  does  not  satisfy  a  variational  principle.  Furthermore, 
LCD  calculations  require  more  computer  time  than  comparable  calculations 
done  using  many-body  perturbation  theory.  However,  CCD  is  more  reliable 
t  an  MBPT  when  more  than  a  single  reference  function  is  important  in  the 
calculation  of  the  electronic  energy.  This  point  has  been  demonstrated 
tor  the  case  of  the  lEg  potential  curve  of  Be2-4 
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APPENDIX  C.  MANY-BODY  PERTURBATION  THEORY 

The  development  of  molecular  many-body  perturbation  theory  (MBPT) 
followed  the  introduction  of  the  occupation  number  representation  ("second 
quantization") ,  from  which  the  diagrammatic  description  of  ordinary  time- 
independent  perturbation  theory  follows. 1  Chapters  five  through  eight 
in  Raimes'  book  outline  the  important  elements  of  the  theory  and  describe 
the  important  linked-diagram  theorem:  only  linked  diagrams  contribute 
to  the  energy  of  the  ground  state. 2  Since  the  literature  of  many-body 
perturbation  theory  is  extensive,  only  an  overview  is  given  here. 3 

Prerequisite  to  developing  MBPT  is  the  partitioning  of  the  mathe¬ 
matical  space  in  which  the  molecular  problem  is  defined.  The  mathematical 
space  for  the  N-electron  problem  is  partitioned  into  the  zero-order  refer¬ 
ence  function,  Do,  and  the  remaining  configurations,  |h)>,  chosen  to  be 
orthogonal  to  D0,  so  that 

P  =  |h><h|h>-1<h| .  (C-l) 

Assuming  the  Moller-Plesset  separation  of  the  molecular  Hamiltonian^  the 
defining  equations  of  perturbation  theory  become. 


H  =  Hq  ♦  v. 

(C-2) 

V  =  E  r^-Zve££(a), 
a<3  ap  a 

(C-3) 

[h(l)  ♦  ve££(l)]  p(l)  -  epXp(l), 

(C-4) 

Hq  =  2[h(i)  +  ve££(i)]. 

(C-5) 

i 


Cl 

S .  Raimes 3  Many-Eleetron  Theory ,  North-Holland^  New  York  (1972). 

C2 

J.  Goldstone,  Proa.  Roy.  Soc.  (London) ,  A229,  267  (1957). 

C3 

R.  J.  Bartlett  and  G.  D.  Purvis ,  "Many-Body  Perturbation  Theory ,  Coupled- 
Pair  Many-Electron  Theory ,  and  the  Importance  of  Quadruple  Excitations 
for  the  Correlation  Problem" ,  Int.  J.  Quantum  Chem.  14s  561  (1978). 

C4 

C.  Moller  and  M.  S.  Plesset ,  "Rote  on  an  Approximation  Treatment  for 
Many-Electron  Systems' Phys.  Rev.  46,  618  (1934). 
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*0  =  Eo  V 

(C-6) 

■  Oolvlv, 

(C-7) 

=  2  e.  , 

(C-8) 

ESCF  Eo  +  E1‘  (c-9) 

All  excited  electronic  configurations,  represented  by  (h) ,  are  eigen¬ 
functions  of  Hq  because  they  are  composed  of  SCF  orbitals  determined  by 
equation  (C-4) .  The  total  energy,  E,  is  given  by 


E  =  E0  +  E1  +  ECORR’  CC-10) 

where  Eqorr  is  the  correlation  energy.  In  Rayleigh-Schrodinger  pertur¬ 
bation  theory  (RSPT),  the  energy  expression  has  the  form, 

4E  =  E  -  E0  -  Ej  *  Ecorr  -  J0^0|V[(E0-  Hq)-1  P(V  -  4E)]m|*o>.  (C-ll) 

Another  expression  for  the  energy  derives  from  the  linked-diagram  expan¬ 
sion, 


AE 


Z  <$>  |v[(E  -  H  ) 

m=oN  o 1  1  v  o  cr 


-1 


(C-12) 


where  the  subscript  "L"  denotes  the  limitation  to  linked  diagrams.  As 
an  instructive  example,  consider  the  second-order  term;  E2, 

E2  =  <^o  I  vl2l><2ll  CE0  -  Hq)  -1 1  h><h [  V |  $Q>*  •  (C-13) 


Since  only  double-excitations  can  mix  accross  the  Hamiltonian  with  a 
Hartree-Fock  reference  function,  the  second-order  energy  expression 
reduces  to 


The  restriction  to  linked-grapha  removes  the  component  from  V\  at 
second-order >  so  that  only  excited  configurations  can  mix  with  that  pro¬ 
duct.  Consequently  the  use  of  as  the  projection  operator  is  jus¬ 
tified. 
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(C-14) 


E2  =  £K<*olvlhK><hKlVl*o>/tEo  -  «$>■ 


E 


Since  k  =1/4  £j.j£ab  D^j ,  representing  the  double-excitations  from 
the  reference  function,  the  second-order  energy  can  be  written  in  terms 
of  two-electron  integrals  as. 


E2  =  1/4  EijEab* l<ijl|ab>|2/(ci 


(C  —  15) 


where  the  double  bar  indicates  an  antisymmetrized  integral.  This  integral 
is  represented  by  the  graph 


At  second-order,  the  use  of  graphical  representations  adorns  with  elegance 
an  otherwise  straightforward  problem.  At  higher-orders,  however,  the 
derivation  of  the  energy  expression  using  the  operator  equation  becomes 
tedious,  so  the  rules  for  directly  drawing  the  linked-diagrams  simplify 
the  analysis.  Bartlett  and  Silver  have  published  a  useful  discourse 
on  rules  for  constructing  the  linked-diagrams.5 

Formally,  the  linked-diagram  expansion  excludes  unlinked-diagrams 
and  this  guarantees  size-extensivity ,  order-by-order.  However,  theoretical 
methods  based  upon  the  linked-diagram  theorem  do  not  satisfy  a  variational 
principal.  As  implemented  by  Bartlett  and  Purvis,  MBPT  provides  a  power¬ 
ful  tool  to  estimate  correlation  energy  at  a  relatively  low  cost.  The 
method  is  limited  to  problems  for  which  $  is  well-approximated  by  a 
single-determinant  reference  function.  Extension  of  the  technique  to 
permit  the  zeroth-order  function  to  be  a  multi-reference  function  have 
been  derived,  however,  so  this  shortcoming  in  the  theory  should  soon  be 
removed .6 


R.  J.  Bartlett  and  D .  M.  Stiver j  nSome  Apects  of  Diagrammatic  Pertur¬ 
bation  Theory ",  Int .  J.  Quantum  Chem .  95>  183  (1975). 

i 

G.  D.  Purvis  and  R.  J.  Bartlett s  private  communication. 
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APPENDIX  D.  UNIFORM  APPROXIMATION  CRITERIA 


Reliable  predictions  of  potential  energy  hypersurfaces  require  com¬ 
putational  accuracy  that  exceeds  the  Hartree-Fock  SCF  level.  While  doubts 
remain  about  the  extent  and  types  of  electron  correlation  effects  that 
must  be  included  to  obtain  accurate  results,  theoretical  methods  that 
meet  certain  criteria  may  provide  more  uniform  approximations  for  energy 
calculations  over  the  entire  potential  hypersurface. 

Two  important  uniform  approximation  criteria  are  "size-extensivity" 
referring  to  the  proper  dependence  of  energy  and  related  properties  on 
the  size  of  the  system* *,!  and  "size-consistency"  and  separability .2,3,4 
Size  extensivity  is  a  characteristic  of  models  based  upon  the  linked- 
diagram  theorem.  While  size  consistency  and  separability  are  often 
equivalent  to  size  extensivity,  size  consistency  may  impose  additional 
requirements  in  some  cases. 

A  method  is  considered  size  consistent  if  the  energy  of  a  system 
made  up  of  two  subsystems,  A  and  B,  infinitely  far  apart  is  equal  to  the 
sum  of  the  energies  of  A  and  B  computed  separately  by  the  same  method. 

For  closed-shell  systems  dissociating  to  closed-shell  fragments,  an  RHF 
(restricted  Hartree-Fock)  reference  function  satisfies  the  size-consistency 
criteria,  and  size  extensivity  (or  the  absence  of  unlinked  diagrams) 
guarantees  size-consistency  for  a  correlated  model  based  upon  that  refer¬ 
ence  function.  On  the  other  hand,  to  satisfy  size  consistency  when  con¬ 
sidering  separation  to  open-shell  fragments  requires  a  UHF  (unrestricted 
Hartree-Fock)  reference  function.  Since  either  an  RHF  or  a  UHF  reference 
function  can  be  used  to  generate  a  linked-diagram,  size-extensive  expan¬ 
sion,  the  condition  of  size  consistency  applied  to  a  correlated  calculation 
based  on  a  single  determinant  reference  function  imposes  the  requirement 
of  a  UHF-type  reference  function  in_  addition  to  the  use  of  a  linked- 
diagram  expansion. 


D1 

R.  J.  Bartlett  and  G.  D.  Purvis ,  ”Many-Body  Perturbation  Theory ,  Coupled- 
Pair  Many  Electron  Theory ,  and  the  Importance  of  Quadruple  Excitations 
for  the  Correlation  Problem” ,  Int .  J.  Quantum  (Them ,  14^  561  (1978). 

D2 

J .  A.  Pople,  J .  S.  Binkley,  and  R.  Seeger,  " Electron  Correlation  Theories 
and  Their  Application  to  the  Study  of  Simple  Reaction  Potential  Surfaces  ”, 
Int.  J.  Quantum  (Them.,  14,  545  (1978). 

D3 

H.  Primas,  in  Modem  Quantum  Chemistry,  edited  by  0.  Sinanoglu ,  Academic 
Press,  NY,  1965. 

D4 

W.  Kutzelnigg,  in  Methods  of  Electronic  Structures  Theory,  edited  by 
H.  F.  Schaefer,  III,  Plenum,  New  York >  1947. 

* 

The  term  ” size-extensivity”  occurs  in  thermodynamics.  An  extensive 
property  is  proportional  to  the  size  of  a  homogeneous  system.  For  a 
homogeneous  system  of  N  units,  energy  predicted  by  a  size-extensive 
method  equals  N  times  the  energy  predicted  for  a  single  unit. 
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A  significant  consequence  of  size  extensivity  occurs  when  studying 
chemical  reactions,  such  as 


A  +  B  -*■  C  +  D. 


The  heats-of-formation,  AHf,  computed  separately  by  a  size  extensive 
method,  combine  to  equal  the  heat-of-reaction,  AHr, 


AHR  =  AHf(C)  +  AHf(D)  "  AHf (A)  -  AHf(B). 


Were  the  heats-of-formation  for  each  species  obtained  by  a  truncated 
Cl  calculation,  this  simple  addition  does  not  equal  the  heat-of-reaction. 
Instead,  to  obtain  a  proper  AH^  requires  supermolecule  calculations  of 
A-B  and  C-D  at  large  internuclear  separation. 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director ,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name: 

Telephone  Number: 


Organization  Address: 


